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ABSTRACT 

""Time-resolved  flowfleld  measurements  for  a  0.5-meter  di¬ 
ameter,  high  work  transonic  turbine  have  been  completed  in  the 
MIT  Turbine  Blowdown  Facility  (TBF).  ~The>  measurements  were 
taken:  to  determine  the  b 1 ade-to-b 1 ade  total  temperature  pro¬ 

file  for  comparison  with  predictions  from  the  Euler  turbine 
equation;  to  determine  the  effect  of  using  time-averaged 
pressures  to  calculate  turbine  performance;  and  to  provide  a 
complete  set  of  time-resolved  turbine  stage  data.  A  prelimi¬ 
nary  ob Ject 1 ve <€f  the  work^ (given  a  6  kHz  blade  passing  fre¬ 
quency)  was  to  deterlmine  the  frequency response  characteris¬ 
tics  of  the  Instrumentation  -employo^(  to  make  the  flowfleld 
measurements.  A  shock  tube  was  built  for  this  purpose. 

Ttre  measurements  were  taken  with  high-frequency  response 
Instrumentation  including  a  dual -hot-wire  aspirating  probe,  a 
’four-way’  angle  probe,  and  two  cobra  head  total  pressure 
probes  Incorporating  silicon  dlaOTagm  pressure  transducers  k 
(Kullte  models  XCQ-093  and  XCQ-062).  The  probes  were  positi¬ 
oned  0.5  axial  chord  lengths  downstream  of  the  turbine  midspan 
and  traversed  c l r cumf er ent 1  a  1 1 y .  The  properties  measured  In¬ 
clude  total  temperature,  total  pressure,  Mach  number,  and  tan¬ 
gential  flow  angle.  Experiments  were  done  at  125X,  100X  and 

80X  corrected  speeds. 

The  aspirating  probe  Is  found  to  have  a  natural  frequency 
of  15.5  kHz  In  the  test  gas  with  a  damping  ratio  of  0.36;  the 
angle  probe  a  char acter I st 1 c  frequency  of  45  kHz  with  a  set¬ 
tling  time  of  18  usee.  Both  results  are  satisfactory  for  ap¬ 
plication  In  the  TBF.  The  measured  total  temperature  profile 
shows  a  peak-to-peak  variation  of  65  C  ( 20X )  and  a  charac¬ 
teristic  frequency  twice  that  of  the  blade  passing  frequency. 
The  peak  total  temperature  Is  out  of  phase  with  the  peak  total 
pressure,  and  suggests  hot  wakes.  The  predicted  temperature 
profile,  made  with  the  Euler  turbine  equation  and  angle  probe 
data,  varies  at  the  blade  passing  frequency  only  and  Is  more 
In  phase  with  the  pressure.  Further,  It  Is  DC  shifted  upward 
20X  from  the  measured  profile.  The  error  associated  with 
using  the  t 1 me-averaged  pressure  from  a  standard  pitot  probe 
( Instead  of  the  actual  time-resolved  signal)  to  calculate  tur¬ 
bine  adiabatic  efficiency  is  shown  to  be  0.25X. 

Thesis  Supervisor:  Dr.  Alan  H.  Epstein 
Title:  Associate  Professor  of  Aeronautics  and  Astronautics 
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speed  channels  can  be  multiplexed  onto  8  of  the  high  speed 
channels.  Programmable  multlspeed  clocks  control  the  data 
sampling  rate  over  the  test  time.  During  the  test,  the  data 
are  stored  In  a  32  million  byte  solid  state  random  access 
memory.  After  the  test,  the  data  are  transmitted  to  a  DEC  PDP 
11/70  for  reduction,  analysis,  and  storage. 

In  summary,  the  TBF  Is  a  versatile,  rigorously  scaled 
vehicle  for  Investigating  the  performance  and  flowfleld 
char acter 1 st 1 cs  of  large  gas  turbines.  The  facility  Is 
particularly  well  suited  for  the  study  of  unsteady  turbine 
flowfleld  phenomena,  especially  as  outlined  In  the  present 
works’  objective  statement.  The  1 nstr ument 1  on  necessary  to 
Investigate  the  turbine  flowfleld  are  now  discussed. 

2.2  Instrumentation  Description,  Performance,  and  Calibration 
The  aspirating  probe,  angle  probe,  silicon  diaphragm 
pressure  transducers,  and  instrumentation  traverser  are 
discussed  In  the  following  sections.  The  physics  of 
operation,  performance  characteristics,  calibration 

procedures,  and  data  reduction  schemes  are  addressed  for  each 
probe.  First,  however,  the  temporal  and  spatial  resolution 
aspects  of  the  Instrumentation  are  dfscussed  relative  to  the 
characteristic  times  and  lengths  associated  with  the  test 
turbine. 

2.2.1  Instrumentation  Spatial  and  Temperoral  Resolution 

Blade  passing  frequency  of  the  test  turbine  In  the  TBF  Is 
slightly  higher  than  6  kHz  when  operated  at  100X  corrected 
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section  Reynolds  number.  An  ar gon-f r eon- 1 2  gas  mixture 
(freon-12  mass  fraction  of  51X)  Is  used  as  the  facility  test 
gas  and  was  chosen  to  accommodate  turbine  scaling.  This  test 
gas  possesses  a  low  gas  constant  and  therefore  a 
correspondingly  low  sonic  speed.  As  such,  the  design  point 
corrected  speed  (tip  Mach  number)  can  be  achieved  with 
relatively  low  mechanical  speeds.  This  In  turn  lowers  the 
minimum  frequency  response  requirements  of  the 
Instrumentation. 

The  design  point  corrected  speed  and  mass  flow  for  this 
turbine  are  achieved  with  Initial  supply  tank  conditions  of 
478  K  temperature  and  4.2  atm  pressure.  The  Initial  rotor 
speed  Is  approximately  100  rev/sec.  Both  the  pressure  and 
temperature  drop  1 sentr op l ca 1 1 y  as  the  supply  tank  'blows 
down*  during  the  test,  but  the  mechanical  speed  of  the  rotor 
Is  controlled  by  an  eddy  current  brake  to  maintain  the  100X 
value  of  corrected  speed.  The  eddy  brake  serves  to  absorb  the 
energy  extracted  from  the  flow  by  the  turbine  during  the  test. 

Instrumentation  access  Is  provided  by  three  13  cm  wide 
windows  spaced  120  degrees  apart  around  the  outer  case  of  the 
turbine.  Each  window  extends  from  upstream  of  the  vane  row  to 
11  cm  downstream  of  the  rotor.  A  schematic  of  the  test 
section  Indicating  the  location  of  the  access  ports  Is  shown 
In  Figure  2-2.  The  turbine  blade  and  vane  dimensions  are  also 
1 nd 1 cated . 

The  TBF  utilizes  a  digital  data  acquisition  system 
designed  for  100  high  speed,  12  bit  channels  with  maximum 
sampling  rates  of  200  kHz  per  channel.  Eight  groups  of  16  low 
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CHAPTER  2 

THE  EXPERIMENT  --  FACILITY,  INSTRUMENTATION,  AND  TECHNIQUE 

A  brief  description  of  the  Turbine  Blowdown  Facility  (TBF)  Is 
given  In  this  chapter.  Facility  setup,  the  test  turbine,  and 
the  data  acquisition  system  are  all  addressed.  The 
Instrumentation  employed  for  the  experiments,  Including  a 
dua 1 -hot-wl re  aspirating  probe,  a  four-pressure-transducer 
angle  probe,  silicon  diaphragm  pressure  transducers,  and  a 
probe  traverser,  are  discussed  in  detail  along  with  their 
respective  physics  of  operation,  performance  characteristics, 
and  calibration  procedures.  Finally,  the  facility  operation 
procedures  and  test  techniques  are  reviewed. 

2.1  The  Facility,  Test  Turbine,  and  Data  Acquisition  System 

The  TBF  Is  a  short  duration  test  facility  presently 
configured  with  a  0.5  m  diameter  high  work,  transonic  turbine 
of  Rolls-Royce  Limited  design  and  manufacture.  Figure  2-1  Is 
a  schematic  of  this  facility.  The  TBF  rigorously  simulates 
actual  turbine  operating  conditions  over  a  test  time  of 
approximately  0.4  seconds.  Design,  construction  and 
operational  aspects  of  the  facility  are  fully  detailed  In  Cl] 
so  only  a  few  relevant  characteristics  are  noted  here. 

The  facility  is  essentially  composed  of  a  large  supply 
tank  separated  from  the  test  turbine  and  dump  tank  by  a  fast 
acting  valve.  The  supply  tank  Is  heated  by  an  oil  Jacket  to 
achieve  the  proper  temperature  ratios  In  the  test  section. 
The  supply  tank  pressure  Is  set  to  achieve  the  desired  test 
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the  literature  with  emphasis  on  the  three-dimensionality  and 
unsteady  nature  of  the  turbine  flowfleld.  Features  of  the 
turbine  flowfleld  described  above,  as  well  as  the  results  of 
the  work  outlined,  will  be  referred  to  again  throughout  this 
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of  the  rotor  airfoils  and  upstream  vanes.  In  this  regard, 
they  describe  two  flow  regimes  cor r espond 1 ng  to  ’maxlmlm 
Interaction*  or  ’minimum  Interaction’  betweem  the  rotor 
leading  edges  and  upstream  vane  wakes.  In  the  case  of  maximum 
Interaction,  the  vane  wakes  and  rotor  wakes  are  said  to  merge 
with  the  consequence  of  leaving  the  blade  passage  uniform 
except  for  two  strong  secondary  flow  vortices.  In  the  minimum 
Interaction  case,  the  vane  wakes  are  said  to  enter  the  blade 
passage,  thereby  Increasing  Its  non-unlformlty  and  masking  to 
some  degree  the  secondary  flow  vortices. 

In  the  paper  above,  a  nond 1  mens  I ona 1  parameter  In  the 
form  of  a  Strouhal  number  Is  used  to  Indicate  the  degree  of 
Interaction  between  the  vane  and  rotor  flowfleld,  and  will  be 
borrowed  for  use  In  Chapter  4.  This  parameter  Is  defined  as 
foil ows : 


S  -  _U  bx  [1.13 

t  Cx 

where  S  Is  a  reduced  frequency  or  Strouhal  number,  U  the  rotor 
speed,  t  the  vane  gap,  bx  the  rotor  axial  chord,  and  Cx  the 
throughflow  velocity.  The  term  t/U  Is  a  time  scale  associated 
with  the  blade  traversal  of  the  upstream  vane  gap.  The  term 
bx/Cx  Is  a  time  scale  associated  with  the  passage  of  fluid 
through  the  rotor.  The  ratio  of  these  terms,  S,  thus  gives  an 
Indication  of  the  number  of  upstream  wakes  In  the  rotor 
passage  at  any  time. 

In  summary,  the  purpose  of  this  section  has  been  to  give 
a  limited  overview  of  recent  turbine  flowfleld  measurements  In 
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near  the  hub.  The  analyses  presented  In  these  papers  can 
essentially  be  extended  to  wake  transport  phenomena  tn 
turb 1 nes . 

Binder  et  al.  [61  employ  a  Laser-Two-Focus  veloclmeter 
to  make  flowfleld  measurements  and  show  that  the  wake  fluid  Is 
highly  turbulent.  Hodson  C71  mounts  hot  wires  on  a  turbine 
rotor  Itself  to  measure  the  time  dependent  velocity  field 
within  the  rotor  passage  and  concludes  that  the  rotor 
unsteadiness  In  dominated  by  the  convection  of  the  wake  fluid 
through  the  rotor .  Both  papers  tend  to  show  that  wake 
generated  unsteadiness,  at  least  at  the  midspan,  can  be 
explained  with  the  wake  cutting  and  migration  theories 
mentioned  above.  Oldfield  and  Doorly  [8]  show  that  In 
addition  to  the  wake,  shock-wave  passing  can  be  a  major  source 
of  flow  unsteadiness  and  may  also  affect  the  rotor  exit 
b 1 ade-to-b 1 ade  profile  via  the  boundary  layer.  Their 
Schl leren  photographs  Indicate  the  upstream  wakes  migrate  to 
the  blade  suction  side  as  they  pass  through  the  blade  row,  as 
predicted  from  the  results  of  C43. 

Finally,  one  of  the  most  complete  sets  of  turbine 
flowfleld  data  Is  given  by  Sharma,  Butler,  Ooslyn,  and  Ortng 
C 9 ] .  They  give  profiles  of  turbulent  energy,  total  pressure, 
static  pressure,  absolute  velocity  magnitude,  and  radial  and 
c  Ir cumf er ent I  a  1  flow  angle  taken  downstream  of  both  the  vane 
and  blade  row.  They  show  that  the  vane  exit  flowfleld  Is 
steady  but  dominated  by  large  Reynolds  stresses  In  the  wake 
regions.  The  rotor  exit  flowfleld.  In  contrast,  was  found  to 
be  highly  unsteady  and  dependent  upon  the  relative  positioning 
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It  Is  clear  from  the  discussion  above  that  the  turbine 
flowfleld  Is  exceedingly  complex.  In  this  thesis,  however,  we 
are  Interested  only  In  mapping  the  b 1 ade-to-b 1 ade  flowfleld 
properties  at  the  rotor  exit  midspan  and,  therefore, 
unsteadiness  associated  with  the  wake  should  be  the  most 
notable  char acter 1 st 1 c  of  the  measurements.  With  this  In 
mind,  some  recent  turbine  and  relevant  compressor  flowfleld 
measurements  In  the  literature  are  now  summarized. 

The  kinematic  analysis  presented  by  Kerrebrock  and 
Mlkolajczak  [41  Is  a  milestone  toward  understanding  the  effect 
of  upstream  wakes  on  the  flow  Inside  blade  passages.  They 
develop  the  concept  of  a  ’slip  velocity’  In  a  transonic 
compressor  to  show  how  the  wake  fluid  does  not  mix  with  the 
invlscld  core  fluid  but,  rather.  Is  transported  relative  to 
the  Invlscld  fluid  as  It  Is  convected  through  the  blade 
passage.  The  result,  which  Is  well  supported  by  experiment. 
Is  a  redistribution  of  stagnation  enthalpy  across  the  blade 
passage  and  Is  characterized  by  an  excess  of  total  temperature 
In  the  stator  wakes.  This  same  phenomena,  as  shown  In  Figure 
1-1,  will  occur  In  a  turbine  except  that  the  wake  will  be 
transported  toward  the  suction  side  of  the  airfoil  rather  than 
the  pressure  side  as  In  the  compressor.  Epstein  and  Ng  C  5  ] 
have  extended  the  analysis  of  C  4 ]  to  Include  quasi-three 
dimensional  effects  of  the  wakes.  In  particular,  they  show 
that  the  relative  flow  angle  Is  not  necessarily  uniform 
through  the  wake  and  this,  coupled  with  the  wakes 
three-dimensionality,  results  In  the  rotor  wake  having  a  lower 
total  temperature  near  the  tip  and  higher  total  temperature 


PAGE  16 


that  Include  full-stage  flowffeld  measurements  and  to  the  best 
of  the  author’s  knowledge,  none  which  offer  high-frequency 
response  temperature  and  pressure  measurements.  So 
ostensibly,  the  temperature  and  pressure  measurements  given 
later  In  this  thesis  are  the  first  such  published. 

As  a  common  theme,  the  published  turbine  data  reveal  a 
flowffeld  dominated  by  secondary  flows  and  unsteadiness,  even 
In  the  rotor  frame  of  reference  [23.  The  variety  of  secondary 
flow  features  serves  as  an  example  of  this  complexity: 
horseshoe  vortices,  formed  by  the  rolling  up  of  the  endwall 
boundary  layer  at  the  leading  edge  of  a  turbine  blade  or  vane; 
passage  vortices,  which  consist  of  a  leg  of  the  horseshoe 
vortex  that  has  migrated  to  the  suction  side  of  the  airfoil  as 
It  Is  convected  through  the  passage;  and  corner  vortices, 
which  rotate  In  the  opposite  sense  of  the  passage  vortex  and 
are  located  always  In  the  endwa 1 1 / suet  1  on  side  corner  of  the 
passage  C  3 1 .  Additionally,  there  are  tip  vortices  In  the 
rotor  due  to  blade  tip  leakage  plus  the  added  Influence  of  a 
centrlflgal  body  force.  Finally,  separation  bubbles  are 
common  on  both  the  vane  and  blade  suction  side  at  the  airfoil 
leading  edge.  Unsteadiness  Is  also  typical  of  the  turbine 
stage  and,  In  fact,  one  can  show  that  for  a  given  particle  In 
Isentoplc  flow  the  flow  must  be  unsteady  for  the  stagnation 
enthalpy  to  change  £23.  Hence,  by  their  very  nature,  turbines 
possess  unsteady  flowflelds.  Summarized,  turbine  unsteadiness 
Is  associated  with  turbulence,  potential  field  Interactions 
between  the  blade  rows,  Inlet  profile  from  a  combustor ,  and 
wakes  generated  by  both  the  vanes  and  blades. 
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Measurements  taken  In  this  facility  form  the  basis  of  the 
effort  described  In  this  thesis. 

The  objective  of  the  work  reported  herein  Is  to  make 
time-resolved  measurements  of  flowfleld  properties  In  a 
full-turbine  stage.  The  measurements  are  taken:  to  determine 
the  b 1 ade-to-b 1 ade  total  temperature  profile  for  comparison 
with  predictions  from  the  Euler  turbine  equation;  to 
determine  the  effect  of  using  time-averaged  pressures  to 
calculate  turbine  performance;  and  to  provide  a  new  set  of 
time-resolved  turbine  stage  data.  The  properties  measured 
Include  total  temperature  In  particular,  along  with  total 
pressure,  static  pressure,  and  radial  and  tangential  flow 
angle.  The  measurements  are  time-resolved  In  the  sense  that 
data  are  taken  at  rates  high  enough  to  distinguish  the  blade- 
to-blade  flowfleld  structure.  To  accomplish  this  objective, 
the  frequency  response  char acter 1 st  I  cs  of  the  1 nstr umentat 1  on 
described  above,  primarily  the  dua 1 -hot-wire  aspirating  probe 
and  angle  probe,  and  their  suitability  for  making  turbine 
flowfleld  measurements  are  Investigated.  The  scope  Is  limited 
to  making  measurements  at  three  operating  points:  125X,  100% , 
and  80%  corrected  speeds  with  their  respective  corrected  mass 
flow  conditions.  The  measurements  are  taken  only  at  the  rotor 
midspan  but  with  circumferential  translation. 

1 . 2  Background 

There  are  an  enormous  number  of  papers  In  the  literature 
describing  both  measured  and  predicted  flowfleld 
characteristics  In  turbine  cascades.  But  there  are  much  fewer 


matched  similarity  parameters,  and  which  utilize  high 
frequency  response  Instrumentation,  are  much  less  common, 
however.  Flowfleld  data  from  the  latter  facilities  are  very 
11ml  ted . 

One  such  facility  which  rigorously  simulates  engine 
conditions  (with  a  simulation  capability  extending  to  40  atm 
Inlet  pressure  and  2500  K  turbine  Inlet  temperature)  has  been 
recently  developed,  constructed,  and  tested  at  the 
Massachusetts  Institute  of  Technology  Cll.  This  facility 
operates  as  a  short  duration  blowdown  tunnel,  having  test 
times  of  approximately  0.4  seconds,  and  Is  configured  with  a 
0.5  meter  diameter,  high  work,  transonic  turbine  designed  and 
manufactured  by  Rolls-Royce  Limited.  High-frequency  response 
Instrumentation  have  also  been  developed  and  complement  the 
utility  of  the  MIT  blowdown  facility.  These  Include  a 
dua 1 -hot-wl re  aspirating  probe  which  yields  a  simultaneous 
measurement  of  flowfleld  total  temperature  and  total  pressure 
In  the  frequency  range  of  DC  to  21  kHz  In  air,  and  an  angle 
probe  Incorporating  four  silicon  diaphragm  pressure 
transducers.  The  angle  probe  provides  measurements  of 
flowfleld  static  pressure,  total  pressure  and  radial  and 
circumferential  flow  angle  and  has  a  frequency  response  of  DC 
to  45  kHz.  Independent  measures  of  total  pressure  can  be  made 
with  cobra  head  Impact  probes.  Together,  the  MIT  blowdown 
facility  coupled  with  the  high-frequency  response 
Instrumentation  provide  a  flexible  vehicle  for  Investigating 
the  complex  flowfleld  of  aircraft  gas  turbines  and  extending 
the  limited  pool  of  turbine  flowfleld  Information. 
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CHAPTER  1 

THE  INTRODUCTION  --  OBJECTIVE,  SCOPE,  AND  BACKGROUND 
This  thesis  deals  with  making  time-resolved  flowfleld 
measurements  In  a  full  gas  turbine  stage.  Properties  measured 
include  total  temperature,  total  pressure,  stattc  pressure, 
Mach  number,  and  tangential  flow  angle.  In  this  Chapter ,  some 
motivation  and  background  for  the  measurements  are  given.  In 
Chapter  2,  the  test  vehicle  employed  for  the  measurements,  the 
MIT  Turbine  Blowdown  Facility,  and  its  high-frequency  response 
instrumentation  are  described.  Chapter  3  gives  the  design  and 
operating  procedures  for  a  shock  tube  facility  constructed  to 
measure  the  frequency  response  characteristics  of  the 
Instrumentation  employed  for  the  turbine  measurements.  In 
Chapter  4,  the  measurements  taken  In  the  turbine  facility  are 
presented  and  briefly  analyzed.  And  finally.  In  Chapter  6, 
conclusions  from  the  work  accomplished  are  summarized. 

1.1  Thesis  Objective  and  Scope 

The  aircraft  high  pressure  turbine  flowfleld  Is 
characteristically  three-dimensional,  unsteady,  and 

compressible,  hence  making  It  exceedingly  difficult  to 
describe  analytically.  Accordingly,  to  refine  performance  or 
emplrlclze  design  codes,  the  turbine  designer  must  rely  on 
experiments  to  determine  the  structure  of  this  spatially 
non-uniform  flowfleld.  Turbine  cascade  data  are  relatively 
abundant  since  facility  costs  are  low  and  the  flowfleld  Is 
benign  compared  to  actual  turbine  operating  environments. 
Facilities  employing  full-stage  test  turbines  with  rigorously 
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w  -  frequency  In  Hertz 

wn  -  undamped  natural  frequency  In  Hertz 
X  -  mole  fraction 

x  -  ellipse  coordinate 

y  -  ellipse  coordinate 

z  -  damp  1 ng  ratio 


Greek 

7^  -  turbine  adiabatic  efficiency 

}j  -  test  gas  viscosity 

p  -  density 

)*■  -  ratio  of  specific  heats 

0  -  abolute  frame  of  reference 

B  -  relative  frame  of  reference 

Super  scr 1 pt 

*  -  reference  variable.  Incompressible  variable 


Subscr 1 pt 
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2  -  second  component  of  binary  gas  mixture,  angle 

probe  transducer  2 
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4  -  turbine  exit  station,  angle  probe  transducer  4 

1  -  refers  to  freestream  value 

n  -  refers  to  local  value 

e  -  tangential  direction 

s  -  refers  to  static  conditions 
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speed.  The  wakes  generated  by  turbine  blades  are  thin,  being 
on  the  order  of  twice  the  blade  trailing  edge  thickness.  In 
fact,  data  obtained  by  Ooslyn,  Drlng  and  Sharma  C  9 ]  suggest 
turbine  blade  wakes  occupy  less  than  10X  of  the  rotor  blade 
spacing.  Accordingly,  to  measure  flowfleld  property 
variations  In  the  region  of  the  wake  (or  In  some  similarly 
small  region  where  property  gradients  are  large)  the 
1  nstr umentat Ion  employed  must  have  frequency  responses  well 
above  the  blade  passing  frequency.  In  the  present  case,  a 
frequency  response  of  something  greater  than  6  kHz  Is  the 
bottom  limit  for  distinguishing  b 1 ade-to-b 1 ade  variations 
while  ten  times  6  kHz  Is  necessary  to  distinguish  property 
variations  through  a  wake.  Hence  the  temporal  resolution  of 
the  Instrumentation  used  Is  a  critical  Item  which  one  must 
reliably  determine.  For  these  reasons,  a  shock  tube  facility 
was  designed  and  constructed  to  Investigate  the  frequency 
response  characteristics  of  the  Instrumentation  used  In  the 
TBF .  The  shock  tube  design  and  operating  procedures  along 
with  test  results  for  the  aspirating  probe  and  angle  probe  are 
presented  (n  Chapter  3. 

Likewise,  the  spatial  resolution  of  the  Instrumentation 
Is  critical  since  their  characteristic  dimensions  are 
essentially  the  dimensions  of  the  smallest  flowfleld  structure 
measureable  by  that  Instrumentation.  And  as  mentioned  above, 
length  scales  associated  with  the  b 1 ade- to-b 1 ade  flowfleld 
structure  are  relatively  small,  especially  In  the  region  of 
the  wake.  We  would  therefore  like  very  small  sensors  and 
probes.  The  size  of  the  aspirating  probe  and  angle  probe 
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relative  to  the  TBF  turbine  blades  and  vanes  are  shown  In 
Figure  2-3.  The  blade  trailing  edge  span  Is  48.4  mm,  the 
blade  spacing  Is  27.9  mm,  and  trailing  edge  thickness  Is 
approximately  0.8  mm.  The  corresponding  character  I st  1  c 
lengths  of  the  aspirating  probe  and  angle  probe  are  3.0  mm  and 
3.3  mm  respectively.  Since  turbine  wake  thicknesses  are  on 
the  order  of  twice  the  trailing  edge  thickness,  we  cannot 
distinguish  Internal  wake  flow  structure  with  the  probes  at 
hand.  We  can  however  distinguish  b 1 ade-to-b 1 ade  structure 
since  both  probes  occupy  only  11X,  or  less,  of  the  blade 
spacing.  From  here  on  we  shall  consider  each  probes  spatial 
resolution  acceptable  for  making  measurements  In  the  TBF. 

2.2.2  The  Dua 1 -Hot-W 1 r e  Aspirating  Probe 

The  aspirating  probe  provides  a  simultaneous  measurement 
of  total  pressure  and  total  temperature  In  an  unsteady 
compressible  flowfleld.  This  probe,  she /n  In  Figure  2-4,  was 
developed  to  Investigate  the  unsteady  flowfleld  downstream  of 
a  transonic  fan  [10]  and  has  the  advantage  of  high  frequency 
response.  Its  development  and  full  description  are  detailed 
In  [11].  Application  of  the  aspirating  probe  to  a  turbine 
environment  Is  a  direct  extension  of  Its  previous  employment 
In  a  transonic  compressor;  however,  certain  aspects  of  the 
probe’s  performance  differ  between  the  two  applications.  In 
particular,  the  probe’s  resolution  and  frequency  response  are 
different.  These  differences  stem  from  the  varying  test  gases 
and  flow  conditions  encountered  with  the  respective  component 
test  facilities,  along  with  choice  of  hot  wire  dimensions. 
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The  aspirating  probe  as  configured  for  application  In  the  TBF 
Is  now  described  along  with  a  review  of  Its  operating 
principles.  Subsequently,  the  probe’s  performance 
characteristics,  calibration  procedures  and  results,  and  data 
reduction  procedures  are  given. 

2.2.2. 1  Description  and  Physics  of  Operation 

A  complete  description  of  the  aspirating  probe  Is  given 
In  [111  but  Is  summarized  here  for  reference.  The  aspirating 
probe  consists  essentially  of  two  coplanar  hot  wires  placed  In 
a  1.5  mm  diameter  constant  cross-sectional  area  channel,  as 
Indicated  In  Figure  2-4.  The  convergent  exit  of  the  channel 
leads  to  a  vacuum  source  and  Is  choked  during  probe  operation. 
The  two  hot  wires  are  set  at  different  overheat  ratios,  and 
hence  different  wire  temperatures,  and  In  this  mode  allow  the 
simultaneous  measurement  of  both  total  temperature  and  total 
pressure,  as  will  be  shown.  A  total  pressure  Impact  probe 
Incorporating  a  Kullte  XCQ-093  silicon  diaphragm  pressure 
transducer  Is  Installed  ’piggyback’  to  the  aspirating  probe. 
The  pressure  transducer  gives  a  redundant  measure  of  total 
pressure  and  allows  for  separate  measures  of  total  temperature 
from  each  hot  wire.  The  pressure  transducer  has  the  added 
benefit  of  reducing  the  probe’s  measurement  uncertainty  but  at 
the  expense  of  reduced  longitudinal  spatial  resolution  till. 
The  hot  wires  themselves  are  made  from  platinum  plated 
tungsten  with  the  sensing  length  defined  by  copper  plating. 
The  wires  are  soldered  In  place  using  multicore  solder  beads. 
Care  Is  taken  to  ensure  the  wires  are  taut  and  as  close  to  the 
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center  of  the  channel  as  possible.  A  typical  wire  separation 
distance  Is  250  urn.  Typical  wire  dimensions  are  1.27  mm  In 
length  and  5.1  urn  In  diameter.  The  hot  wires  are  operated 
with  standard  TSI-1050  constant  temperature  anemometers. 

The  classical  governing  equation  for  the  anemometer 
bridge  voltage  of  a  constant  temperature  hot  wire  Is: 

V2  -  Frs  +  Rwl2  Jl  1  k  Nu  [2.1] 

tTw  -  r  T1  Rw 

where  V  Is  the  bridge  output  voltage,  Tw  the  hot  wire 
temperature,  r  the  wire  recovery  factor,  Rs  the  anemometer 
resistance  In  series  with  the  hot  wire,  Rw  the  wire  operating 
resistance,  1  the  wire  length,  k  the  fluid  thermal 
conductivity,  and  Nu  the  Nusselt  number.  This  equation 
applies  only  when  the  hot  wire  Is  operated  In  a  gas  of  uniform 
composition.  The  Nusselt  number  Is  a  function  of  many  fluid 
and  wire  properties  C121,  but  Is  often  expressed  In  the 
following  form  as  a  function  of  the  hot  wire  Reynolds  Number 
on  1  y  i 

Nu  -  a  (Raj"  ♦  b  [2.2] 


The  coefficients  of  this  equation,  which  Is  often  referred  to 
as  King's  Law,  are  determined  by  calibration.  Since  the 
channel  of  the  aspirating  probe  Is  choked,  we  may  represent 
the  channel  mass  flux  as  a  function  of  the  fluid  properties, 
total  pressure,  total  temperature,  and  ratio  of 


son  1 c-to-channe 1  area  ratio  at  the  hot  wire  plane: 


*+1 


The  corresponding  Reynolds  number  Is  then  given  by: 


[2.3] 


[2.4] 


where  d  Is  the  wire  diameter,  the  fluid  viscosity,  P  the 
total  pressure,  T  the  total  temperature,  A*  the  choked  exit 
area,  Ac  the  wire  plane  channel  area,  y  the  fluid  ratio  of 
specific  heats,  and  R  the  fluid  gas  constant.  Finally,  we  may 
combine  Equations  2.1,  2.2,  and  2.4  to  yield  the  governing 
equation  for  the  aspirating  probe: 


2 

,  J'+'l 

N  \ 

V2 

=.  f  Rs+Rw]  n  1  k 

)a  d[  P  4: 

*  b 

[2.5] 

[Tw  -  r 

T]  Rw 

)  ]T[-fP  Ac  \|  R 

s 

which  can 

be  reduced  to  the 

V2 

following  form. 

c  «0"  * D 

[2.6] 

[Tw  -  r  T] 

Equat 1  on 

2 . 6  app lies  when 

the  probe  geometry 

and 

gas 

composition  are  fixed.  By  operating  the  two  wires  at 


different  overheat  ratios  In  separate  anemometer  circuits,  two 
forms  of  Equation  2.6  are  obtained  from  which  we  may  determine 
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the  unknown  total  pressure  and  total  temperature.  The 
coefficients  C  and  D  are  found  by  calibration  and, 
unf or tunate 1 y ,  are  a  function  of  gas  temperature.  This 
functional  relationship  Is  driven  by  the  sensitivity  of  the 
test  gas  viscosity  and  thermal  conductivity  to  temperature. 
This  characteristic  of  the  aspirating  probe  Is  discussed 
further  In  Section  2. 2. 2. 3. 

The  Influence  of  the  hot  wire  properties  on  the  probe 
output  voltage  Is  evidenced  by  Equation  2.5.  In  particular, 
increasing  wire  resistance,  length,  and  dtameter  all  tend  to 
Increase  the  anemometer  output  voltage.  The  net  effect  Is 
generally  Improved  resolution.  This  aspect  of  the  aspirating 
probe  along  with  Its  other  performance  character  1 st 1 cs  are  now 
cons  1 der ed . 

2.2.2. 2  Performance  Character  1 st  1  cs 

As  Indicated  above,  the  hot  wire  characteristics 
including  length,  diameter,  and  resistance  dictate  the 
aspirating  probe  output  voltage  for  a  given  gas.  Nominal  wire 
conditions  of  1.27  mm  length,  5.1  um  diameter,  and  3.00  Ohms 
resistance  give  a  temperature  sensitivity  of  12.5  mv/C  at  0.86 
atm  and  an  overheat  ratio  of  1.8.  At  an  overheat  ratio  of 
2.0,  the  temperature  sensitivity  Is  10.1  mv/C  at  0.86  atm.  At 
a  temperature  of  350  K,  the  probe’s  sensitivity  to  pressure  Is 
0.88  V/atm  at  an  overheat  ratio  of  1.8  and  1.00  V/atm  at  an 
overheat  ratio  of  2.0.  These  senslvltles  apply  to  an 
ar gon-f r eon- 1 2  mixture  typical  of  that  used  In  the  TBF  and, 
further,  the  pressure  of  0.86  atm  and  temperature  of  350  K  are 
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character 1 st 1 c  of  the  flowfleld  downstream  of  the  facility's 
turbine. 

The  hot  wires  are  operated  with  two  separate  TSI-1050 
constant  temperature  anemometers.  The  noise  levels  are 
typically  on  the  order  of  2  mv  peak-to-peak .  The 
ana  1 og-to-d 1 g  i  ta 1  converter  employed  has  a  resolution  of  2.5 
mv,  and  as  such,  sets  the  probe  resolution  at  approximately 
0.2  C  in  temperature  and  0.003  atm  in  pressure.  An  error 
analysis  at  1  atm  and  290  K  gives  a  pressure  uncertainty  of 
0.4X  and  a  temperature  uncertainty  of  0.4X  (1.20  when  data 
are  reduced  with  1  wire  and  the  pressure  transducer  C201.  For 
operation  in  dynamic  flowflelds,  angle  effects  are  negligible 
below  15  degrees  and  the  response  is  symmetric  about  the  probe 
axis  til].  Finally,  for  the  TBF  test  gas,  the  aspirating 
probe  natural  frequency  is  15.46  kHz  with  a  damping  ratio  of 
0.36.  The  determination  of  frequency  response  is  detailed  in 
Chapter  3. 

2. 2. 2. 3  Calibration  Procedures  and  Results 

The  aspirating  probe  Is  calibrated  against  temperature 
and  pressure  In  an  argon-f reon- 1 2  gas  mixture  with  a  freon 
mass  fraction  of  approximately  0.51  X.  The  temperatures  and 
pressures  surveyed  Include  the  operating  range  of  the  test 
turbine.  The  calibrations  are  accomplished  in  the  TBF  supply 
tank.  Ideally,  we  would  like  to  calibrate  the  probe  while 
immersed  In  a  steady  flow  (or  well  characterized  dynamic  flow) 
but  no  facility  for  accomplishing  this  type  of  calibration  was 
available.  The  static  calibration  Is  fully  acceptable. 
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though,  given  the  high  frequency  response  of  the  probe. 

The  coefficients  C  and  D  are  determined  from  a 
least-squares  linear  regression  of  Equation  2.6,  resubmitted 
below,  for  a  set  of  voltage  and  pressure  pairs  at  a  constant 
calibration  temperature  and  hot  wire  overheat  ratio. 


T3 


C2.6] 


The  exponent,  N,  Is  chosen  and  kept  constant  for  the  curve 
fit.  Iteration  gives  a  value  of  N,  usually  around  0.3, 
yielding  the  lowest  standard  deviation  versus  the  calibration 
data.  Typical  values  for  the  aspirating  probe  coefficients 
are  given  for  two  hot  wire  overheat  ratios  In  Table  2-1  below. 
Hot  wire  properties  and  dimensions  are  also  Indicated. 


Overheat  Ratio 

1.80 

2.00 

C  CT**(N/2-l )  V 

**2  /  P**N3 

0.3115 

0.2752 

D 

[ V**2  /  T] 

-0.0231 

-0.0226 

N 

0.30 

0.30 

Calibration  Gas 

Argon  Freon 

-12  ml xture 

with 

Freon  mass 

fraction  of 

0 . 5 1 X 

W 1  re  Mater  1  a  1 

Platinum  coated  Tungsten 

Wire  Temperature, 

Tw  C  K  3 

491 . 1 

539.4 

Wire  Length,  L 

[MM3 

1  .27 

1.27 

Wire  0 lameter ,  D 

C  uM3 

5.  10 

5.  10 

Temperature  Coefficient  of 

Resistance,  a 

C0HMS/C3 

0.004 

0.004 

T  -  Temperature 

In  Kelvin 

P  -  Pressure  In 

Atmospheres 

V  -  Volts 

Table  2-1*  Aspirating  Probe  Calibration  Coefficients 
Versus  Hot  Wire  Overheat  Ratio 


PAGE  32 


The  gas  mixture  kinematic  viscosity,  |4 •  and  thermal 
conduct l v l ty ,  k,  are  both  functions  of  temperature  and  thus 
have  the  unfortunate  effect  of  making  C  and  D  functions  of 
temperature  as  well.  The  functional  relationships  between  C 
and  0  and  the  fluid  properltes  are  clear  from  Equation  2.5. 
From  Collls  and  Williams  [13],  the  fluid  properties  take  on 
values  commensurate  with  the  average  of  the  freestream  gas 
temperature  and  the  hot  wire  operating  temperature.  This 
explains  the  variation  of  coefficients  C  and  D  with  overheat 
ratio,  and  hence  wire  temperature.  In  Table  2-1 .  As  long  as 
the  overheat  ratio  Is  kept  constant,  though,  the  wire 
temperature  effect  Is  not  a  problem.  The  coefficient 
sensitivity  to  the  freestream  gas  temperature  on  the  other 
hand  Is  a  problem  since  It  Is  this  temperature  that  we  wish  to 
measure.  Accordingly,  the  coefficients  must  be  corrected  for 
fluid  property  variations  with  temperature. 

Figures  2-5  A  and  B  show  calibration  curves  for  the 
aspirating  probe  at  overheat  ratios  of  1.8  and  2.0 
respectively.  Two  sets  of  calibration  data  are  shown  on  each 
figure  for  gas  freestream  temperatures  of  321.7  K  and  349.5  K. 
The  data  taken  at  321.7  K,  Indicated  by  the  squares,  are  curve 
fitted  to  Equation  2-6  yielding  the  coefficients  shown.  For 
purposes  of  Illustration,  these  same  coefficients  along  with 
Equation  2.9  are  used  to  predict  the  second  set  of  calibration 
data  taken  at  349.5  K,  Indicated  by  the  triangles.  For  both 
overheat  ratios  the  standard  deviation  Increases  from 
approximately  0.3  C  In  the  curve  fitted  case  to  4.0  C  In  the 
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predicted  case. 

Presumably,  we  can  Improve  the  prediction  by  accounting 
for  the  effects  of  freestream  temperature  on  the  gas  mixture 
viscosity  and  thermal  conductivity.  Such  corrections  would 
take  the  following  form  where  the  primed  variables  Indicate 
reference  values. 


C 

D 


o  *  _k_ 

k » 


C2.71 

C  2 . 8  ] 


Note 

we  may 

use  this  same  approach  to 

correct  for 

changes 

1  n 

the 

wire 

characteristics  like 

res  1  stance , 

1 ength , 

and 

d 1 ameter  . 

These  corrections  are 

straightforward  and 

not 

addressed  further . 

As  previously  noted,  the  fluid  properties  are  evaluated 
at  the  average  of  the  freestream  temperature  and  hot  wire 
operating  temperature.  Values  for  the  kinematic  viscosity  and 
thermal  conductivity  are  taken  from  C14J  for  argon  and  CIS] 
for  freon-12.  The  empirical  correlation  proposed  by 
Buddenberg  and  Wilke  C16],  developed  to  predict  the  constants 
In  the  Suther 1 and-Th 1 esen  equation  for  binary  gas  mixtures.  Is 
used  to  determine  the  properties  of  the  ar gon-f reon- 1 2  test 
gas  mixture.  This  correlation  Is  shown  below  for  the  test  gas 
mixture’s  thermal  conductivity;  the  same  correlation  applies 
for  kinematic  viscosity. 
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k  mix 


kl  ♦  k  2 


1  ♦  X2  Q12 

5TT 


1  ♦  XI  Q21 

71 


12.3  ] 


where 


Q12 


Q21 


[2. 10] 


[2.11] 


In  the  equations  above,  1  and  2  refer  to  each  component  of  the 
binary  mixture  with  X  representing  mole  fraction  and  M 
molecular  weight. 

Applying  the  corrections  above  results  In  a  much  Improved 
prediction  of  the  calibration  data  as  Is  shown  In  Figures  2-5 
C  and  D.  These  data  are  the  same  as  discussed  previously.  In 
each  case,  the  standard  deviation  of  the  predicted  data  Is  of 
the  order  of  the  curve  fitted  data.  Clearly,  the  method  for 
correcting  the  aspirating  probe  calibration  coefficients 
works.  Hence,  using  some  Iteration  scheme  to  find  C  and  D  as 
functions  of  temperature  relative  to  the  calibration  data,  we 
may  solve  for  both  total  temperature  and  total  pressure  In  the 
turbine  flowfleld.  The  aspirating  probe  data  reduction 
procedures,  Including  this  Iteration  scheme,  are  addressed  In 
the  following  section. 
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2. 2. 2. 4  Data  Reduction  Procedures 

The  three  analog  signals  from  the  two  hot  wires  and  one 
total  pressure  probe  are  digitized  during  the  test  and 
subsequently  stored  on  hard  disk  In  a  POP  11/70  computer.  The 
silicon  diaphragm  pressure  transducers  provide  good  measures 
of  total  pressure  so  we  are  primarily  Interested  In  measuring 
total  temperature  with  the  aspirating  probe.  There  are  two 
approaches  to  reducing  the  data. 

The  first  approach  utilizes  the  total  pressure 
measurement  made  with  the  ’piggybacked’  pressure  transducer. 
The  pressures  measured  by  the  attached  probe  and  voltage  from 
either  wire  Is  used  to  calculate  temperature  from  Equation 
2-6.  For  this  calculation,  the  exponent  N  Is  chosen  based  on 
experience  with  the  calibration  data.  A  simple  secant-type 
Iteration  scheme  Is  then  employed  to  solve  for  temperature 
while  Including  the  effects  of  the  temperature  on  the 
coefficients  C  and  D  C173.  The  value  of  tungsten’s 
temperature  coefficient  of  resistance,  e,  Is  Input  to 
calculate  the  wire  temperature  with  the  equation  belowi 


Rh  ■  Overheat  Ratio 
Rc 

=  1  +  e  C  Tw  -  r  T 1  [2.12] 

where  Rh  Is  the  hot  wire  operating  resistance,  Rc  the  hot  wire 
cold  resistance,  Tw  the  wire  temperature,  r  the  fluid  recovery 
factor,  and  T  the  fluid  temperature.  Both  the  wire  overheat 
ratio  and  recovery  factor  are  Input.  For  thin  wires  the 
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recovery  factor,  r,  Is  always  greater  than  0.98  C121. 

The  second  approach  utilizes  both  wires  to  find 
temperature;  the  pressure  transducer  Is  not  used.  This 
approach  takes  advantage  of  the  Coefficient  correction  scheme 
addressed  In  the  previous  section.  Based  on  the  calibration 
data,  curves  of  bridge  output  voltage  for  both  wires  can  be 
constructed  versus  temperature  and  pressure  as  shown  In  Figure 
2-6.  Using  the  numerical  equivalent  of  this  grid  and  a  simple 
Interpolation  routine,  we  can  estimate  the  flowfleld 
temperature  given  the  voltage  signals  from  the  two  hot  wires. 

2.2.3  The  Four-Pressure-Transducer  Angle  Probe 

The  angle  probe  Is  a  high-frequency  response  Instrument 
which  provides  measurements  of  total  pressure,  static 
pressure,  Mach  number,  and  radial  and  tangential  flow  angles, 
and  Is  shown  In  Figure  2-7.  As  with  the  aspirating  probe, 
this  probe  was  designed  and  constructed  to  Investigate  the 
flowfleld  downstream  of  a  transonic  compressor  C101.  The 
original  angle  probe  was  spherical  In  shape  and  Included  five 
silicon  pressure  transducers  from  which  the  flowfleld  data 
were  retrieved  C18H191.  Subsequent  designs  Incorporated  four 
pressure  transducers.  Instead  of  five,  mounted  on  an 
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2.2.3. 1  Description  and  Physics  of  Operation 

The  angle  probe  consists  of  four  silicon  diaphragm 
pressure  transducers  located  on  an  elliptical  stem  as  shown  In 
Figure  2-7.  The  transducers  utilized  are  Kullte  model  XCQ-093 
with  abolute  pressure  reference,  and  are  numbered  on  the 
figure  for  reference.  Transducers  2  and  3  are  offset  from  the 
center  transducer  by  45  degrees  while  transducer  4  Is  placed 
flat  on  the  probe  tip.  The  tip,  as  shown,  Is  also  set  at  a  45 
degree  angle.  Cooling  water  fs  supplied  through  the  base  of 
the  probe  and  circulated  around  the  backside  of  the 
transducers,  thus  enhancing  their  thermal  stability. 

As  was  done  with  the  spherical  probe  [19],  we  may  define 
coefficients  based  upon  the  four  pressures  measured  by  the 
angle  probe.  These  ares 
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and  results  given.  Finally,  measurements  taken  In  the  TBF 
coupled  with  their  brief  analysis  are  presented. 
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traverse  velocity  by  the  300  ms  point.  The  full  traverse 
takes  place  over  approximately  120  ms. 

All  the  probes  described  above  operate  throughout  the 
full  test  time  with  the  exception  of  the  aspirating  probe. 
This  probe’s  hot  wires  handle  dynamic  gas  loads  well  but  not 
particulates  In  the  test  gas.  The  wires  are  susceptible  to 
particulate  deposits  which  may  weaken  the  wires  or  otherwise 
Invalidate  their  calibrations.  Though  the  facility  test  gas 
Is  considered  clean,  care  Is  taken  to  operate  the  probe’s 
suction  only  for  the  test  time  (note  the  wires  are  left  on 
continuously).  A  small  pneumatically  activated  fast  acting 
valve  Is  used  for  this  purpose.  The  valve  Is  spliced  Into  the 
flexible  hose  connecting  the  aspirating  probe  to  a  small 
mechanical  vacuum  pump.  The  valve  Is  opened  at  150  ms  and 
flow  through  the  probe  Is  Immediately  established.  The  valve 
closes  100  ms  later  but  It  takes  approximately  4  seconds  for 
the  flow  Into  the  aspirating  probe  to  stop.  The  probe  does 
not  unchoke  until  the  1  second  point,  well  after  the  test 
time.  By  this  method,  dangers  to  the  aspirating  probe’s  hot 
wires  are  minimized. 

In  summary,  we  have  reviewed  the  TBF  setup  and  operating 
procedures,  paying  particular  attention  to  the  operation  of 
the  aspirating  probe.  Elsewhere  in  this  chapter,  the  TBF 
facility,  turbine,  and  data  acquisition  system  have  been 
described.  The  various  probes  and  probe  traverser,  along  with 
their  descriptions  and  modes  of  operation,  have  also  been 
addressed.  In  the  following  chapters,  a  shock  tube  facility 
for  measuring  Instrumentation  frequency  response  Is  described 
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of!  through  the  flow  control  valve  and  against  the  drive 
piston.  The  magnitude  and  steadiness  of  the  piston  velocity 
Is  controlled  by  both  the  gas  cylinder  pressure  and  control 
valve  setting.  The  probes  are  attached  directly  to  the  drive 
piston.  The  travel  of  the  device  Is  60  mm  and  takes  place 
nominally  In  0.1  seconds.  A  typical  traverse  Is  depicted  in 
Figure  4- 1 . 

2.3  Facility  Operation  Procedures  and  Experimental  Technique 

Operation  of  the  TBF  Is  now  reviewed.  A  more  detailed 
discussion  of  the  facilities  operating  procedures  are  given  In 
Cl].  Once  the  Instrumentation  required  for  a  test  is 
Installed,  the  valve  Is  closed  between  the  supply  and  dump 
tanks  and  both  are  evacuated  (<0.1  torr).  The  supply  tank  Is 
then  heated  to  the  desired  temperature,  typically  478  K.  The 
heating  Is  done  slowly  to  ensure  thermal  equilibrium.  Once  at 
temperature,  the  supply  tank  Is  filled  with  the  test  gas 
mixture  of  argon-freon-12  to  the  desired  pressure,  typically 
4.2  atm.  The  gas  assumes  the  supply  tank  temperature  by  the 
end  of  the  fill.  Finally,  the  rotor,  which  Is  still  In 
vacuum.  Is  brought  up  to  speed  by  a  small  electric  motor. 
When  the  desired  Initial  rotor  speed  Is  achieved,  the  valve 
and  eddy  current  brake  are  simultaneously  triggered  and  flow 
of  the  test  gas  from  the  supply  tank  Is  established. 
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transducers  which  are  normally  referenced  to  vacuum. 
Calibrations  are  performed  Immediately  before  and  after  each 
test  by  alternating  the  reference  pressure  from  vacuum  to 
atmospher  e . 

2.2.5  The  Traverser  Mechanism 

Short  test  times  char acter  1  s  1 1 1  c  of  the  TBF  complicate 
the  radial  and  c 1r cumf erent f a  1  traversal  of  tnstrumentat (on 
during  a  test.  The  traversal  must  be  timed  to  occur  during 
the  high-speed  data  acquisition  phase  and,  preferably,  take 
place  with  constant  velocity  to  facilitate  data  reduction. 
Cao  [241  designed  and  constructed  such  a  traverser  for  use  In 
the  TBF;  It  Is  detailed  In  Figure  2-13.  The  traverser  can 
move  probes  both  circumferentially  and  radially  downstream  of 
either  the  rotor  or  vane  row.  Circumferential  traverses  alone 
were  utilized  In  the  present  set  of  experiments  so  only  this 
function  of  the  traverser  Is  described. 

The  c 1 r cumf erent (a  1  traversing  mechanism  has  as  Its 
essential  component  a  hydraulically  driven  piston.  The 
device’s  high  pressure  cylinder  Is  charged  with  argon  (or 
other  gas)  to  approximately  21  atmospheres.  This  cylinder,  by 
way  of  a  solenoid  valve.  Is  connected  to  an  oil  cylinder. 
This  In  turn  (s  connected  to  the  device’s  drive  piston  through 
a  flow  control  valve  (used  to  maintain  constant  linear 
velocity.  A  linear  differential  variable  transformer  (LDVT) 
attached  to  the  piston  Is  used  to  Indicate  position. 

Actuation  of  the  device  occurs  when  the  solenoid  valve  is 
triggered,  thus  allowing  the  high  pressure  gas  to  force  the 
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The  scheme  then  corrects  theta  and  phi  for  a  new  Mach  number 
calculated  from  the  new  static  pressure  and  total  pressure 
calculated  from  Equation  2.25.  The  coefficients  H23,  CPI,  and 
KP2  (or  KP3)  are  subsequently  updated  and  the  scheme  continues 
until  convergence  of  the  static  pressure  Is  achieved.  Upon 
convergence  the  Mach  number  Is  re-calculated  based  upon  the 
test  gas  ratio  of  specific  heats  as  previously  Indicated.  The 
final  values  of  the  current  Iteration  are  then  stored  and  used 
as  Initial  guesses  for  the  follow-on  data  set.  Typically, 
five  Iterations  are  required  before  convergence  Is  achieved. 

2.2.4  The  Silicon  Diaphragm  Total  Pressure  Transducers 

Two  silicon  diaphragm  pressure  transducers,  Kullte  models 
XCQ-093  and  XCQ-062,  provide  high-frequency  total  pressure 
measurements.  The  XCQ-093  Is  carried  piggyback  with  the 
aspirating  probe  and  Is  similar  to  the  ones  used  In  the  angle 
probe.  With  a  diameter  of  approximately  1.6  mm  (0.032 
Inches),  the  XCQ-062  subminiature  pressure  transducer  Is  the 
smallest  Kullte  pressure  transducer  available  providing  both 
static  and  dynamic  flowfleld  measurements.  The  XCQ-062  has  a 
natural  frequency  In  excess  of  500  kHz  and  has  an  operating 
temperature  range  of  -55  C  to  120  C.  Its  nominal  sensitivity 
Is  26.5  mv/atm.  Both  transducers  are  constructed  similarly  to 
pitot  probes  with  cobra  head  designs,  as  shown  In  Figure  2-4 
for  the  XCQ-093.  Each  transducer  diaphragm  Is  recessed 
approximately  one  diameter  from  the  cobra  head  entrance  to 
negate  the  probe's  sensitivity  to  flow  angle  variations  In  the 
♦20  degree  range.  Both  total  pressure  probes  use  differential 
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the  tangential  angle,  theta,  Is  found  by  linear  Interpolation. 
Initial  guesses  for  the  total  and  static  pressure  are  made 
with  PI  taken  as  the  total  pressure  and  the  lesser  of  P2  and 
P3  as  the  static  pressure.  From  these  an  'air  equivalent’ 
Mach  number  Is  calculated  using  the  Isentroplc  pressure  ratio 
relation.  The  term  air  equivalent  Is  used  to  Indicate  that 
the  Mach  number  Is  based  on  air’s  value  of  the  ratio  of 
specific  heats  since  the  calibrations  were  done  In  air.  Upon 
convergence  of  the  Iteration  scheme,  the  true  value  of  the 
Mach  number  based  upon  the  test  gas  properties  Is  found.  The 
next  coefficient  calculated  Is  CP4  which  in  turn  Is  used  to 
find  phi  as  a  function  of  theta.  Having  Initial  values  for 
Mach  number,  theta,  and  phi  the  coefficients  H23,  CPI,  and  KP2 
(or  KP3  depending  on  which  of  P2  or  P3  was  used  as  an  Initial 
guess  for  static  pressure)  are  determined  by  Interpolation. 
The  static  and  total  pressures  are  then  calculated  with  the 
following  relations  based  upon  the  coefficient  definitions: 


Ps  =  P 2 , 3  -  P2P3  KP 2 , 3  12.241 
Pt  =>  Ps  +  H  <  1  -CP  1  >  P2P3  12.251 
where  P2P3  *  (P2-P1)  ♦  (P3-P1)  C2.261 

If  the  static  pressure  calculated  from  Equation  2.24  differs 
from  the  fnltal  assumed  static  pressure,  a  new  value  of  static 
pressure  Is  chosen.  This  new  value  of  static  pressure  Is 
taken  as  the  average  of  the  assumed  and  calculated  pressures. 
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test  flowfield.  Coefficient  F23  is  least  sensitive  to  the 
radial  angle  showing  virtually  no  effect  until  theta  reaches  .+ 
20  degrees.  The  remaining  coefficients  show  consistent  but 
only  moderate  reactions  to  phi. 

In  summary,  the  calibration  results  are  supported  by  the 
simple  prediction  model  detailed  above.  It  is  clear,  however, 
that  the  angle  probe’s  unique  geometry  influences  the 
calibration  coefficients.  It  is  noted  that  the  pneumatic 
model  of  the  angle  probe  used  for  the  calibrations  is 
suppossed  to  be  an  accirate  replica  of  the  instrument  used  In 
the  TBF  .  If  variations  exist,  the  calibration  data  may  not  be 
fully  applicable.  Unfortunately,  comparative  calibrations 
were  not  accomplished.  The  geometry  related  aspects  of  the 
angle  probe  are  Important  and,  accordingly,  are  Included  in 
the  data  reduction  schemes.  These  are  now  addressed. 

2. 2. 3. 4  Data  Reduction  Procedures 

The  angle  probe  yields  a  rather  complex  set  of 
coefficients  Influenced  both  by  the  flowfield  and  probe 
geometry,  as  was  seen  in  the  calibration  section.  The  data 
reduction  scheme  outlined  below  has  been  used  for  both  the 
original  spherical  probe  and  earlier  f our-tr ansducer  models. 
It  relies  heavily  upon  linear  interpolation  and,  thus,  all  of 
the  geometry  related  nuances  of  the  calibration  data  are 
retained.  The  procedure  for  the  current  probe  and  most  of  the 
data  reduction  software  has  been  developed  by  Gertz  1231,  and 
is  summarized  In  Figure  2-12. 

The  scheme  begins  with  the  calculation  of  F23  from  which 
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since  In  each  case  the  coefficients  vary  conslstantly  with  the 
progression  of  the  two  Independent  parameters.  The  peculiar 
’bumps*  and  shifts  In  the  curves  are  taken  to  be  associated 
with  the  unique  geometry  of  the  angle  probe  and  are  accounted 
for  In  the  data  reduction  schemes. 

Angle  probe  coefficients  F23,  H23,  KP2,  CPI,  and  CP4  are 
depicted  In  Figures  2-10A  through  2-10E  versus  tangential 
angle  for  Mach  numbers  of  0.35,  0.5,  0.6,  and  0.85.  The 
calibration  data  were  taken  for  zero  radial  angle. 
Coefficient  F23  Is  relatively  unaffected  by  Mach  number  except 
In  the  +8  degree  regions.  The  ’bumps’  at  these  locations  are 
not  born  out  by  the  flowfleld  prediction  and  are  most  likely 
tied  Into  the  probe  geometry.  Coefficient  H  Is  strongly 
affected  by  Mach  number  via  the  total  pressure,  as  expected. 
Coefficient  KP2  (as  well  as  KP3)  Is  also  strongly  affected 
and,  as  seen  with  F23,  shows  a  peculiar  anomaly  at  -8  degrees 
which  becomes  more  pronounced  with  decreasing  Mach  number. 
Again,  this  Is  most  likely  associated  with  the  placement  of 
the  pressure  transducer  on  the  probe  shaft  and  Its  attendant 
flat  face.  KP3,  though  not  shown,  exhibits  the  same  behavior 
as  KP 2 . 

The  same  angle  probe  coefficients  as  above  are  depicted 
In  Figures  2-11A  through  2-11E  versus  tangential  angle.  For 
these  coefficients,  the  Mach  number  Is  set  at  0.6  and  the 
radial  angle  Is  varied  from  -10  degrees  to  ♦10  degrees  In  5 
degree  Increments.  Coefficient  CP4  Is  strongly  affected  by 
the  radial  angle  variation  and,  as  discussed  In  the  data 
reduction  section,  Is  key  to  retrieving  this  property  of  the 
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to  the  expected  turbine  midspan  exit  Mach  number.  The 
prediction  and  data  agree  well  considering  the  difference  In 
the  actual  probe  geometry  and  that  of  an  Ideal  two-dimensional 
elliptical  cylinder.  In  particular,  the  transducer  faces  on 
the  angle  probe  are  flat  and  thus  maintain  a  higher  stagnation 
pressure  than  would  an  Ideal  ellipse  at  the  same  location. 
This  Is  evidenced  by  CPI  where  the  calibration  data  Is 
conststantly  higher  than  the  prediction.  This  coefficient’s 
Insensitivity  to  theta  values  In  the  _+  5  degree  range  Is 
clear.  Likewise,  the  data  and  predictions  for  CP2  and  CP3  are 
of  the  same  form  as  CPI  but  offset  by  the  transducer 
separation  angles  of  +.45  degrees  from  center.  Coefficient  F23 
Is  determined  from  PI,  P2,  and  P3  In  such  a  way  that  the 
difference  In  the  calibration  data  and  prediction  Is 
nullified,  resulting  In  a  suprlslngly  good  fit.  Coefflcents 
H23  and  KP2,  unlike  F23,  Include  the  freestream  static 
pressure  In  their  definitions  and  therefore  show  a  difference 
beteen  the  measured  and  predicted  cases.  For  both 
coefficients  the  data  are  lower  than  the  predictions  but 
clearly  show  Identical  trends.  The  difference  can  be  traced 
to  the  larger  values  of  PI,  P2  and  P3  used  to  normalize  the 
coefficients.  KP3,  the  mirror  Image  of  KP2,  exhibits  the  same 
behavior  as  KP2.  Based  on  these  results,  the  method  for 
calibrating  the  probe  In  the  freejet  Is  deemed  acceptable. 
This  being  the  case,  the  effects  of  Mach  number  and  radial 
angle  on  the  angle  probe  coefficients  are  now  shown.  These 
calibration  data  are  presented  to  Illustrate  the  effects  of 
Mach  number  and  radial  angle;  only  brief  remarks  are  made 
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these  velocities,  we  can  solve  for  the  local  Mach  number,  Mn, 
given  the  freestream  Mach  number,  Mt,  as  followst 


Mn*  -  Un2  [2.221 

{t  R  Tn) 


Finally,  we  may  solve  for  the  local  stat 1 c-to-tota 1  pressure 
ratio  using  the  standard  tsentroptc  form  as  a  function  of 
local  Mach  number i 

i 


Ln 

Pt 


■  ♦  v- 1  Mn2J 


[2.23] 


Thus  we  can  determine  the  local  static  pressure  at  each 
transducer  location  and,  except  for  CP4,  predict  all  the  angle 
probe  coefficients  once  the  freestream  Mach  number  and 
tangential  flow  angle,  theta,  are  specified  and  the  radial 
angle,  phi.  Is  set  at  zero. 

Measured  and  predicted  angle  probe  coefficients  F23,  H23, 
KP2,  CPI,  CP2,  and  CP3  for  a  Mach  number  of  0.6  and  zero 
radial  angle  are  shown  In  Figures  2-9A  through  2-9F.  This 
Mach  number  was  chosen  for  the  predictions  since  It  Is  close 
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range  of  calibration  Mach  numbers  Is  0.35  to  0.85  so 
compressibility  effects  should  be  considered.  For  this  case 
we  can  apply  a  subsonic  Prandtl -Glauert  correction  as  follows 
where  the  primed  variables  apply  to  the  Incompressible  case 
1221  » 


B 

X 


[2.20] 


Applying  this  correction  to  Equation  2.19  and  rearranging 
yields i 


[2.21  ] 


Thus,  Equation  2.21  provides  the  local  tangential  velocity 
distribution  over  an  elliptical  cylinder  once  given  the 
ellipse  fineness  ratio,  A/B,  and  the  freestream  velocity.  We 
are  only  Interested  In  the  local  velocity  at  each  transducer 
relative  to  the  freestream  angle  of  attack  and,  therefore, 
require  definitions  for  determining  the  appropriate  Internal 
angles,  Q.  These  definitions,  shown  In  Figure  2.8,  apply  only 
to  a  circular  cylinder  but  the  error  In  applying  the 
prediction  to  the  angle  probe  Is  small  since  It  has  a  low 
fineness  ratio  (A/B  Is  appr ox Imatel y  1.35).  Additionally,  the 
results  only  apply  when  the  the  radial  angle  Is  zero.  Using 
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blockage  effects  were  small.  Given  this  result,  the  small 
nozzle  was  utilized  for  all  the  angle  probe  calibrations  and 
allowed  measurements  to  be  recorded  over  the  full  Mach  number 
range  of  Interest. 

To  Increase  confidence  In  the  calibration  data, 
predictions  can  be  made  for  the  coefficients  F23,  H23,  KP3, 
KP2,  CPI,  CP2,  and  CP3  by  modelling  the  angle  probe  as  a 
two-dimensional  elliptical  cylinder.  By  using  the  potential 
solution  for  such  an  elliptical  cylinder,  we  can  calculate  the 
local  surface  velocity  and,  given  the  freestream  total 
temperature,  determine  the  local  Mach  Number.  Finally,  we  can 
determine  the  ratio  of  local  stat 1 c-to-tota 1  pressure  ratio  at 
each  transducer  location,  thus  allowing  direct  calculation  of 
the  coefficients.  The  equations  for  this  procedure  now 
follow. 

The  potential  solution  for  the  surface  flow  velocity  over 
a  two-dimensional  elliptical  cylinder  Is  [211: 
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The  axes  and  coordinate  system  are  defined  In  Figure  2-8.  The 


nominal  thermal  zero  shifts  are  on  the  order  of  3X  per  full 
scale  per  55.5  C  1100  F)  while  the  nominal  sensitivity  shifts 
are  less  than  1.5%  per  full  scale  per  55.5  C  (100  F>. 
Finally,  the  compensated  temperature  range  Is  275  K  to  422  K 
( 35  to  300  F ) . 

2. 2. 3. 3  Calibration  Procedures  and  Results 

The  angle  probe  Is  calibrated  In  a  freejet  against  radial 
angle,  tangential  angle,  and  Mach  number.  The  calibration 
test  gas  Is  air.  The  calibrations  are  accomplished  with  a 
precisely  duplicated  pneumatic  model  of  the  Instrument  used  In 
the  TBF .  This  Is  done  to  avoid  damaging  the  relatively 
fragile  and  expensive  probe  constructed  with  the  silicon 
diaphragm  pressure  transducers.  All  of  the  calibrations  are 
done  statically.  As  discussed  later  in  Chapter  3,  we  would 
like  to  know  the  response  of  the  angle  probe  to  dynamic 
flowflelds  but  no  means  are  currently  available  to  produce 
we  1 1 -char acter I  zed  dynamic  flowflelds  with  frequencies  above  5 
kHz  [201.  Accordingly,  only  the  static  response  of  the  angle 
probe  has  been  Incorporated  Into  the  data  reduction  schemes. 

Probe  blockage  errors  must  be  considered  when  freejets 
are  used  for  calibrations.  Accordingly,  two  freejet  nozzles 
of  2.54  cm  and  5.08  cm  diameter  were  employed  for  the  angle 
probe  calibrations.  Due  to  facility  constraints,  the  maximum 
steady  Mach  number  achievable  with  the  large  nozzle  was  only 
0.4,  well  below  the  desired  maximum  calibration  value  of  0.85. 
Comparston  of  coefficient  values  taken  by  both  nozzles  at  the 
same  Mach  number  showed  no  differences,  suggesting  that  probe 
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Pt  and  Ps  refer  to  the  freestream  total  and  static  pressures 
respectively.  Two  additional  coefflcents,  CP2  and  CP3,  can  be 
defined  similarly  to  CPI.  These  coefficients  are  not  used  In 
the  data  reduction  but  are  considered  In  the  calibration 
section  presented  later. 

Based  on  the  definitions,  F23  and  CP4  provide  directional 
sensitivity,  H23  Mach  number  sensitivity,  and  KP2  and  KP3 
static  pressure  sensitivity.  Coefficient  CPI  aids  In 
determining  the  total  pressure.  In  theory  then,  we  may  use 
the  set  of  coefficients  to  retrieve  flowfleld  property 
Information  Including  the  radial  and  tangential  flow  angles, 
total  pressure,  static  pressure,  and  Mach  number  once  the 
probe  Is  fully  calibrated. 

2. 2. 3. 2  Performance  Characteristics 

We  are  particularly  Interested  In  the  frequency  response 
of  the  angle  probe.  The  Kullte  model  XCQ-093  pressure 
transducers  used  In  the  Instrument  have  natural  frequencies  In 
excess  of  200  kHz.  However  In  practice,  the  length  scale 
characteristic  of  the  probe,  Its  diameter,  coupled  with  the 
flowfleld  velocity  sets  the  characteristic  frequency  of  the 
probe.  This  frequency  Is  approximately  45  kHz  and  therefore 
well  above  blade  passing  frequency.  The  angle  probe’s 
frequency  response  characteristics  are  addressed  further  In 
Chapter  3. 

The  pressure  transducers  are  approximately  2.4  mm  (0.093 


Inches)  In  diameter  and  have  nominal  sensitivities  of  31.6 
mv/atm  with  a  nonlinearity  of  less  than  0.5%  full  scale.  The 


CHAPTER  3 


THE  SHOCK  TUBE  --  MEASUREMENT  OF  PROBE  FREQUENCY  RESPONSE 

This  chapter  details  the  development  and  construction  of 
a  shock  tube  test  facility  built  for  the  purpose  of  measuring 
Instrumentation  frequency  response.  Predictions  and 
measurements  of  frequency  response  are  given  for  both  the 
aspirating  probe  and  angle  probe. 

3.1  Shock  Tube  Description  and  Operating  Procedures 

The  design  approach  and  description  of  the  shock  tube 
facility  are  detailed  In  this  section.  The  operating 
procedures,  shakedown  results,  and  data  acquisition  system  are 
also  presented. 

3.1.1  Shock  Tube  Design  Approach  and  Desrlptlon 

A  shock  tube  test  facility  was  designed  and  constructed 
for  the  purpose  of  measuring  Instrumentation  frequency 
response.  In  particular,  the  frequency  response 
characteristics  of  the  aspirating  probe  and  angle  probe  are 
desired  since  they  both  operate  downstream  of  the  TBF  turbine 
rotor.  In  this  turbine  environment,  as  mentioned  In  Chapter 
2,  the  probes  are  subjected  to  high  frequency  flowfleld 
fluctuations  which  occur  mostly  at  the  blade  passing  frequency 
and  Its  harmonics.  In  order  to  accurately  measure  or  even 
distinguish  the  b 1 ade-to-b 1 ade  flowfleld  structure,  the  probes 
natural  frequencies  must  be  well  above  the  blade  passing 
frequency.  Therefore,  we  must  determine  each  probes  frequency 


response  In  order  to  assess  Its  applicability  for  use  In  the 
turbine  facility.  The  shock  tube  Is  a  vehicle  for  making  this 
assessment . 

Step  function  tests  are  commonly  used  to  Investigate 
Instrumentation  dynamic  response.  As  a  general  rule  for  step 
testing,  the  rise  time  of  the  step  function  must  be  less  than 
one-fourth  of  the  natural  period  of  the  Instrument  tested 
[253.  This  rule  ensures  excitation  of  the  Instruments  natural 
oscillations.  Shock  tubes  typically  provide  very  sharp 
changes  (on  the  order  of  10  A~8  seconds)  In  pressure, 
temperature,  and  flow  velocity,  thus  making  them  suitable  for 
Investigating  natural  frequencies  well  above  one  megahertz. 
Further,  they  can  be  designed  for  simple  operation  while  still 
providing  repeatable  results.  It  was  for  these  reasons  that  a 
shock  tube  test  facility  was  selected  for  the  frequency 
response  measurements. 

In  approaching  the  shock  tube  design  task  several 
criteria  were  considered,  the  first  being  ease  of  operation. 
In  particular  It  was  felt  the  shock  tube  should  operate 
without  a  diaphragm  bursting  mechanism.  This  was  a  point  of 
concern  since  the  anticipated  test  pressure  ratios  were  low. 
Accordingly,  any  variance  In  the  diaphragm  material  properties 
could  result  In  very  different  burst  pressures.  Also, 
bursting  could  be  ragged,  thus  sending  particles  of  diaphragm 
material  Into  the  driven  tube  which  could  potentially  damage 
sensors.  Secondly,  the  facility  had  to  be  flexible  In  the 
sense  that  any  practical  bottled  gas  could  be  used  In  either 
the  driver  or  driven  sections  of  the  tube.  This  criterion 
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required  that  fill  and  vacuum  lines  be  sent  to  both  sides  of 
the  shock  tube.  Next,  the  device  had  to  provide  repeatable 
pressure  rises  on  the  order  of  the  b 1 ade- to-b 1 ade  pressure 
fluctuations  downstream  of  the  turbine  rotor.  This  criterion 
essentially  dictated  the  diaphragm  pressure  ratio.  Finally, 
the  test  time  --  the  time  beginning  with  the  passage  of  the 
Initial  shock  wave  and  ending  with  the  passage  of  either  the 
reflected  shock  or  gas  Interface  --  had  to  be  long  enough  to 
allow  for  several  periods  of  a  low  order  kilohertz  signal. 
This  criterion  dealt  primarily  with  the  driven  tube  length  and 
relative  positioning  of  the  Instrumentation  ports. 

The  shock  tube  was  designed  In  accordance  with  the 
criteria  above.  Both  the  driver  and  driven  sections  of  the 
device  were  set  at  lengths  of  3.281  m  (10  ft),  and  were 
constructed  from  shedule  40  (2  In  l.d.)  pipe  made  of  316 
stainless  steel.  The  opposite  ends  of  each  section  were 
capped  with  304  stainless  steel  blind  flanges  while  the  mating 
ends  were  capped  with  similar  socket  welded  flanges.  A  double 
O-rlng  and  groove  configuration  was  machined  Into  the  mating 
flanges.  This  arrangement  proved  successful  In  holding  the 
diaphragms  taut  while  the  driver  side  of  the  facility  was 
charged.  No  diaphragm  breaking  mechanism  was  required,  as 
discussed  later.  Add  1 1 Iona  1 1 y ,  stainless  steel  dowel  pins 
were  placed  Into  the  mating  flanges  to  help  with  alignment. 
The  pins  also  served  to  hold  the  diaphragms  In  place  while  the 
facility  was  recycled  between  tests.  To  allow  for  operation 
w.th  any  bottled  gas,  fill  and  vacuum  lines  were  lead  to  the 
opposite  ends  of  each  section.  Standard  one-quarter  Inch 


copper  tubing  was  used  for  the  fill  lines;  one-half  Inch 
copper  tubing  for  the  vacuum  lines.  Between  the  driver  and 
driven  sections,  the  fill  lines  were  connected  with  a  flexible 
length  of  PVC  tubing  while  the  vacuum  lines  were  connected 
with  flexible  vacuum  hose.  Valving  was  designed  to  allow 
filling  and  evacuation  of  either  section  Independently  from 
the  other.  Finally,  facility  Instrumentation  was  minimized 
and  Included  only  a  Type  0  thermocouple  and  vacuum  and 
pressure  gauges. 

Positioning  of  the  Instrumentation  mounts  was 
accomplished  by  predicting  the  shock  velocity,  reflected  shock 
velocity,  and  Interface  gas  velocity  as  a  function  of 
diaphragm  pressure  ratio,  gas  temperature,  and  gas  properties. 
Standard  one-dimensional  shock  tube  relations  were  employed 
for  these  calculations  C263.  A  map  of  time  versus  driven  tube 
position  can  be  constructed  by  dividing  the  velocities  above 
by  the  chosen  driven  tube  length,  as  shown  In  Figure  3-J  for 
air.  The  diaphragm  pressure  ratio  for  the  case  shown  Is 
approximately  1.48,  yielding  a  pressure  ratio  across  the  shock 
of  approximately  1.2.  This  pressure  rise  Is  on  the  order  of 
what  an  instrument  would  experience  downstream  of  the  TBF’s 
turbine  rotor.  The  Initial  values  for  temperature  and 
pressure  In  the  driven  tube  are  standard  room  conditions.  As 
previously  defined,  the  test  time  begins  with  the  passing  of 
the  Initial  shock  and  Its  attendant  step  change  In  pressure, 
temperature,  and  velocity.  The  test  time  concludes  with  the 
passage  of  either  the  reflected  shock  wave  or  Interface  gas 
velocity  and  Its  respective  muddling  of  the  original  step 
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function  flow  field.  For  the  conditions  of  Figure  3-1,  the 
maximum  test  time  occurs  at  the  24X  driven  tube  length  and  has 
a  duration  of  approximately  13  milliseconds.  This  Is  an 
enormous  amount  of  time  for  dynamic  testing  of 
Instrumentation.  Higher  diaphragm  pressure  ratios  would  give 
higher  shock  velocities  and,  for  the  same  driven  tube  length, 
have  the  dual  effect  of  reducing  overall  test  time  and  moving 
the  maximum  test  time  position  farther  down  the  tube.  In 
anticipation  of  using  higher  pressure  ratios  at  some  future 
date,  the  Instrumentation  mount  locations  were  set  at  the  30X 
and  40X  driven  tube  positions.  Test  times  for  these  positions 
In  air  are  12  and  10  milliseconds  respectively. 

The  next  step  following  the  design  and  construction  of 
the  shock  tube  Involved  the  determination  of  diaphragm 
materials  and  operating  procedures.  These  aspects  of  the 
shock  tube  development  along  with  Initial  shakedown  results 
are  now  presented. 

3.1.2  Shock  Tube  Operating  Procedure  and  Shakedown  Results 

Operation  of  the  shock  tube  proved  to  be  very  simple. 
For  tests  In  air,  the  driven  tube  Is  simply  filled  with  air 
from  a  regulated  high  pressure  cylinder  until  the  diaphragm 
breaks.  Depending  on  the  diaphragm,  this  process  was 
surprisingly  repeatable.  For  operation  with  other  gases,  both 
sides  of  the  shock  tube  are  Initially  evacuated  and  then 
filled  with  the  test  gas  to  some  pre-determ  1 ned  pressure.  The 
driven  section  fill  line  Is  then  cut  off  from  the  supply  while 
the  driver  continues  to  fill  until  the  diaphragm  breaks.  This 
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process  also  proved  to  be  repeatable.  Separate  gases  could  be 
used  In  the  driver  and  driven  sections  but  this  permutation 
was  not  Investigated.  In  all  cases,  the  double  O-rlng  and 
groove  configuration  of  the  mating  flanges  held  the  diaphragms 
secure  until  bursting  occurred. 

Several  potential  diaphragm  materials  were  tested 
Including  cellophane,  paper,  vellum,  and  metal  foil. 
Diaphragms  that  broke  repeatably,  cleanly,  and  at  a  pressure 
ratio  of  1.2  to  2.0  (normalized  by  atmosphere)  were  desired. 
Cellophane  (DuPont  MSD-60,  220  gauge,  22.8  pm  thick)  was  found 
to  meet  these  desirable  qualities  the  best.  This  material  did 
have  the  occasional  problem  of  releasing  particles  upon 
breaking,  but  no  sensors  were  damaged.  The  repeatability  on 
breaking  was  approximately  15%  (typically  7  psla  _+  1  psla, 
which  was  within  the  accuracy  of  the  pressure  gauge  used). 

Shock  wave  velocity  measurements  were  conducted  to 
confirm  the  device  operated  as  standard  shock  tube  relations 
would  predict.  The  wave  velocities  were  measured  using  two 
platinum  thln-fllm  gauges  operated  with  TSI-1050  constant 
temperature  anemometers.  The  thln-fllm  gauges  were  placed  on 
the  shock  tube  Inner  wall  surface  at  a  known  separation 
distance.  The  time  between  response  to  a  passing  shock  wave 
was  measured  on  an  oscilloscope  and  used  along  with  the 
separation  distance  to  calculate  the  wave  velocity.  The 
velocities  measured  were  consistently  within  10%  of 
predictions.  These  results  were  considered  good  given  that 
the  diaphragm  break  pressure  was  read  visually  from  a  pressure 
gauge  mounted  to  the  driver  tube.  These  results  were  later 
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confirmed  using  a  hot  wire  to  measure  the  post-shock  flow 
field  velocity  and  Its  characteristic  steadiness  over  the  test 
time.  These  data  are  discussed  later  In  Section  3.3.2. 
Additionally,  It  was  found  the  wave  velocities  from  the 
diaphragm  pressure  ratio  desired,  approximately  1.5,  were 
sufficient  to  excite  the  damped  natural  frequencies  of  Kullte 
XCQ-093  and  XCQ-062  strain  gauge  pressure  transducers.  These 
frequencies  were  250  kHz  and  500  kHz  respectively,  well  above 
those  anticipated  from  either  the  aspirating  probe  or  angle 
probe . 

A  Tekronlcs  468  digital  recording  oscilloscope  Is  used  to 
acquire  the  test  data.  The  measuring  Instrument  acts  as  the 
trigger.  Once  taken,  the  data  are  passed  to  an  LSI  11/23 
computer  through  an  IEEE-488  GPIB  and  stored  on  floppy  disk. 
The  data  are  subsequently  transferred  to  a  DEC  PDP  11/70 
computer  for  reduction  and  analysis. 

In  summary,  a  shock  tube  facility  for  testing 
Instrumentation  frequency  response  was  designed  and 
constructed.  Its  operation  was  very  simple  while  Its  results 
proved  to  be  repeatable.  Measured  shock  wave  velocities  and 
Interface  gas  velocities  were  consistently  within  10X  of 
predictions  from  standard  one-dimensional  shock  wave 
relations.  With  the  shock  tube  operating  as  a  viable  tool  for 
frequency  response  Investigations,  measurements  were  made  with 
both  the  aspirating  probe  and  angle  probe.  These  results  are 
now  presented. 

3.2  Aspirating  Probe  Frequency  Response 
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In  this  section,  the  governing  equations  for  the 
aspirating  probe  natural  frequency  and  damping  ratio  are 
described.  Subsequently,  we  employ  the  equations  to  predict 
the  probe’s  natural  frequency  and  damping  In  both  air  and 
Freon-12,  and  compare  the  predictions  to  shock  tube  test 
results.  The  Instrument’s  predicted  response  In  an 
argon-freon  mixture,  the  gas  utilized  In  the  TBF ,  Is  then 
presented.  Finally,  we  discuss  means  of  Improving  the 
frequency  response  characteristics  of  the  aspirating  probe 
based  on  the  shock  tube  results. 

3.2.1  Governing  Equations  for  Natural  Frequency  and  Damping 

The  frequency  response  of  compensated  constant 
temperature  hot  wire  anemometers  Is  typically  on  the  order  of 
100  kHz  C121.  In  the  present  case,  however,  frequencies 
associated  with  the  gas  flowing  In  the  aspirating  probe 
channel  must  be  considered  and,  as  will  be  shown  In  the  next 
sections,  are  considerably  lower  than  the  hot  wire  anemometer 
frequency  quoted  above. 

The  aspirating  probe  channel  can  be  modelled  as  a 
constant  area  duct  with  mass  flow  as  shown  In  Figure  3-2.  As 
with  a  number  of  fluid  systems  having  tubing  or  ducts 
(pressure  transducers  with  plumbing,  for  example),  the 
aspirating  probe  can  be  shown  to  exhibit  second-order 
Instrument  behavior.  In  particular,  Whitehead  1271  gives  the 
rel at lonsh Ips  required  to  Investigate  a  system  as  modelled 
above.  These  relationships  are  addressed  In  the  next  section, 
but  first,  a  brief  discussion  of  second-order  Instrument 
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behavior  Is  given. 

Using  the  definitions  of  £281,  a  second-order  Instrument 
Is  defined  as  one  that  follows  the  equation: 

al  dj  +  a2  +  a0  q  *  b0q1  £3.1] 

where  q  Is  the  output  quantity,  qt  Is  the  Input  quantity,  t  Is 
time,  and  the  a's  and  b’s  are  system  physical  parameters  which 
we  shall  assume  to  be  constant.  From  equation  3.1,  we  may 
define  three  parameters: 

K  »  b0  static  sensitivity 

a V 

wn  a  I a0  undamped  natural  frequency, 

^a2  rad/tlme  <  f  ■  w/2  Hertz  > 

z  =  a  1  damping  ratio,  dimensionless 
2  ^a0aZ 

Rearranging  equation  3.1  In  terms  of  these  three  parameters 
yields: 

I'D2  ♦  2zD  +1  1  q0  ■  K  ql  £3.23 

L"wn2  wn  J 

where  D  »  d  /dt. 

From  equation  3.2,  then,  we  may  define  an  operational  transfer 
function  applicable  to  all  second-order  Instruments: 

g0  (D)  -  K _  £3.33 

q  1  Dz  / wn 2  *  2zD/wn  *  1 
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It  Is  evident  from  equation  3.3,  that  to  determine  the 
operating  character s  1 1 1  c s  of  any  second-order  Instrument,  we 
need  only  determine  Its  natural  frequency  and  damping  ratio. 

To  Investigate  the  response  of  a  second-order  Instrument 
to  dynamic  Input  conditions  --  Its  frequency  response  --  we 
may  use  equation  3.3  along  with  a  suitable  sinusoidal  transfer 
function  representing  some  harmonic  excitation.  Following  the 
analysis  In  1291,  for  the  complex  sinusoidal  transfer  function 
q0/q1  ( 1w) ,  we  obtain  from  equation  3.3t 

30  (  iw)  =  K _  [3.4] 

qT  (Iw/w ny*  2z1w/wn  +  1 

where  1  *  -l 

w  =»  frequency, 

which  may  be  reduced  to  the  following  formi 

q0  <1w>  «■  _ 1 _  L  <t>  [3.5] 

Kq  1  £  1  -  (w/wn)2-J*+4zzwi/wn 

where  ■  arctan  2z  [3.6] 

w/wn  -  wn/w 

For  any  chosen  Input  frequency,  w,  the  magnitude  of  equation 
3.5  represents  numerically  the  ratio  of  the  output  to  Input 
signal  amplitudes.  The  angle  of  this  complex  number 
represents  the  phase  angle  by  which  the  output  leads  the  Input 
signal.  The  term  q0/q1K  Is  often  referred  to  as  the 
magnification  ratio  while  the  term  w/wn  Is  called  the 
frequency  ratio. 

Using  equations  3.5  and  3.6,  we  may  construct 
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nond 1 mens  1 ona 1 1  zed  frequency  response  curves  by  choosing 
values  of  frequency  ratio  and  damping  ratio,  as  shown  In 
Figure  3-3.  Curves  of  phase  angle  versus  frequency  ratio  are 
also  shown.  These  curves  represent  the  frequency  response 
characteristics  of  all  second-order  Instruments.  From  the 
curves,  we  see  that  Increasing  natural  frequency,  wn,  yields  a 
greater  range  of  frequencies,  w,  for  which  the  magnification 
ratio  Is  close  to  unity.  Accordingly,  a  high  natural 
frequency  Is  required  to  accurately  measure  high  frequency 
Input  signals.  The  effect  of  damping  ratio,  z,  Is  also 
apparent  as  It  can  be  seen  that  a  value  of  approximately  0.6 
to  0.7  gives  the  widest  flat  magnification  ratio.  In 
addition,  for  these  values  of  damping  ratio,  we  see  that  the 
phase  angle  varies  almost  linearly  with  frequency  ratio.  This 
linearity  Is  Important  If  we  wish  to  obtain  the  correct  shape 
of  the  Input  signal  C291. 

In  summary,  all  second-order  Instruments  are  governed  by 
the  operational  transfer  function  given  by  equation  3.3,  which 
has  as  Its  essential  variables  the  natural  frequency,  wn ,  and 
damptng  ratio,  z.  To  determine  the  frequency  response  of 
these  Instruments,  we  may  use  a  sinusoidal  transfer  function 
which  Is  obtained  by  replacing  Iw  for  D  In  equation  3.3,  thus 
obtaining  equation  3.4.  From  equation  3.4,  we  can  construct 
plots  of  magnification  ratio  and  phase  angle  versus  frequency 
ratio  with  the  damping  ratio  as  an  Independent  variable.  From 
these  plots  It  can  be  seen  that  a  high  natural  frequency  Is 
required  to  accurately  measure  high  frequency  Input  signals. 
Additionally,  It  Is  seen  that  damping  ratios  on  the  order  of 
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0.6  to  0.7  give  the  wfdest  range  of  flat  magn f f f cat f on  ratio 
and  phase  angle  linearity.  Finally,  since  the  aspirating 
probe  behaves  as  a  second-order  Instrument,  we  can  Investigate 
Its  frequency  response  by  determining  Its  natural  frequency 
and  damping  ratio.  This  Is  done  In  the  following  section. 

3.2.2  Prediction  of  Probe  Natural  Frequency  and  Damping  Ratio 
As  mentioned  above,  the  aspirating  probe  can  be  modelled 
as  a  constant  area  duct  with  mass  flow.  For  the  conditions  of 
subsonic  mean  flow  and  a  choked  exit,  Whitehead  [27]  obtains 
the  following  relation  for  the  system  undamped  natural 
frequency  In  Hertz,  f,  as  a  function  of  the  mean  flow  Mach 
number,  M,  gas  sonic  speed,  C,  and  duct  length,  L» 

2 

f  »  (2n-l  >(  I-M  )C  C 3 . 7 ] 

- “TC - 


The  variable  n  Is  an  Integer  value  defining  the  frequency  mode 
shape.  For  the  analysis,  Whitehead  assumes  the  duct  has  a 
uniform  cross  section  and  Is  fed  at  constant  total  pressure 
through  a  short  contraction.  He  also  makes  the  assumptions  of 
one-d Imens Iona  1  flow,  perfect  gas,  small  flow  oscillations, 
and  negligible  fluid  friction.  To  the  extent  that  these 
assumptions  are  realistic  for  the  probe  while  operating 
downstream  of  a  turbine  rotor,  the  result  applies.  The  only 
marginal  assumption  seems  to  be  the  one  regarding  constant 
total  pressure.  In  fact,  total  pressure  will  vary  as  a 
minimum  at  the  blade  passtng  frequency.  Nevertheless,  as  will 
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be  shown,  the  probe’s  natural  frequency  fs  well  above  the 
blade  passing  frequency,  so  on  the  time  scale  of  the  natural 
frequency  the  total  pressure  can  be  considered  quas 1 -steady . 
We  shall  assume  from  now  on  that  Whitehead’s  result  Is 
applicable  and  represents  the  governing  equation  from  which  we 
may  predict  the  aspirating  probe’s  natural  frequency. 

In  addition  to  the  natural  frequency,  Whitehead  gives  the 
following  relationship  for  determining  the  system  damping 
ratio,  z,  as  a  function  of  Mach  number,  M,  and  ratio  of 
specific  heats.  Gamma: 


<2n-l)  z  3  In  (1+MHl+kM)  C3.83 

-  mTTTPETT) 

where  k  =  (Gamma-1) 

- 2 - 

Having  relationships  for  both  natural  frequency  and  damping 
ratio,  we  can  now  make  predictions. 

It  Is  reasonable  but  not  certain  that  equations  3.7  and 
3.8  apply  to  the  aspirating  probe;  but  by  testing  the  probe 
in  two  different  gases,  thereby  changing  the  gas  constant  and 
gamma  and  hence  sonic  speed,  we  can  assess  their 
applicability.  If  It  Is  clear  they  do  apply,  then  we  can 
predict  the  probe’s  frequency  response  for  the  ar gon-f r eon- 1 2 
mixture  used  In  the  TBF  with  confidence.  Measurements  could 
be  made  In  the  appropriate  ar gon-f r eon- 1 2  mixture,  but 
loglstlcally  and  operationally  It  Is  more  convenient  (and 
equally  valid)  to  make  the  measurements  In  two  more  readily 
available  gases.  Accordingly,  air  and  freon-12  were  chosen 
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for  the  shock  tube  tests.  These  gases  have  the  added 
advantage  of  possessing  greatly  different  gas  constants,  thus 
giving  greatly  different  predictions  for  natural  frequency. 
The  predictions  for  these  gases  at  room  temperature  and  for 
the  actual  aspirating  probe  channel  length  are  given  below  In 
Table  3-1. 


Test  Gas 

Air 

Freon- 1 2 

Channel  length,  L 

[MM] 

3.886 

3.886 

Temperature,  T 

CK] 

305.  J0T 

305.0 

Gas  constant,  R 

tJ/KG  K] 

287.2 

68.8 

Gamma 

1  .  40 

1.14 

Sonic  speed,  C 

[M/SJ 

350.2 

146.3 

Channel  Mach  number,  M 

0.  199 

0.2 

Natural  Frequency, 

f  CkHzl 

21.63 

9.53 

Damp  1 ng  Ratio,  z 

0.97 

0.89 

Table  3-1:  Aspirating  Probe  Natural  Frequency 
and  Damping  Ratio  Prediction  for 
Both  Air  and  Freon-12 


The  channel  Mach  number  noted  In  the  table  above  Is  calculated 
using  the  channel  cross-section  to  sonic  orifice  area  ratio 
and  the  ratio  of  specific  heats,  gamma.  In  this  calulatlon, 
growth  of  the  thermal  and  momentum  boundary  layers  are  also 
considered.  The  predictions  can  now  be  compared  to 
measurements  made  for  the  aspirating  probe  In  the  same  set  of 
gases . 


3.2.3  Measurement  of  Probe  Natural  Frequency  and  Damping  Ratio 
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of  course,  Is  due  to  the  probe  translation  across  the  vane 
gap.  Comparison  of  the  7  and  37  cycle  ensembles  Indicate  only 
small  dlsrepancles  despite  the  probe  translation. 
Accordingly,  the  data  presented  In  the  following  sections  are 
ensembled  over  37  cycles  except  where  noted  otherwise. 

4.2.2  Aspirating  Probe  Measurements 

Total  temperature  and  pressure  ratios  measured  with  the 
aspirating  probe  are  presented  In  this  section.  The 
measurements  were  taken  at  the  1 00%,  125%,  and  80X  corrected 
speed  points.  Table  4-3,  below,  gives  a  summary  of  the  test 
conditions  and  reference  values.  Figures  showing  the  test 
data  at  each  test  point  are  also  Indicated. 


Test 

Corrected 

Speed 

T  1  me 

Tref 

£K3 

Pref 

Eatml 

F  1  gur e 

2 

NC  =  100X 

292 

461 

3.43 

4-7,  4-8 

360 

459 

3.36 

4-9,  4-10 

430 

456 

3.26 

4-11,  4-12 

3 

NC  *  125% 

292 

476 

3.50 

4-13 

355 

473 

3 . 37 

4-14 

430 

470 

3.26 

4-15 

4 

NC  *  80X 

292 

468 

3.44 

4-16 

365 

465 

3 . 36 

4-17 

4  30 

462 

3.23 

4-18 

Table  4  —  3  s  Summary  of  Aspirating  Probe  Measurement 
Cond  1 1 1  ons 


The  data  taken  with  the  aspirating  probe  hot  wires  were 
reduced  with  the  pressure  signal  obtained  from  the  XCQ-092 
pressure  transducer,  thus  yielding  two  separate  measures  of 
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value  for  reduced  frequency  Indicates  there  Is  at  least  one 
vane  wake  In  a  blade  passage  at  any  Instant  In  time.  Further, 
this  Implies  that  we  can  determine  the  predominant  rotor  exit 
profile  with  relatively  few  blade-passing  cycles.  This  point 
Is  confirmed  shortly. 

Ensemble  averages  of  total  pressure  are  presented  In 
Figures  4-3  A  and  8  for  the  292  msec  test  point.  Each  period 
of  the  ensembled  traces  shown  are  the  same;  seven  periods  are 
displayed  only  for  comparison  with  the  Instantaneous  traces. 
The  measurements  were  taken  at  the  100X  corrected  speed 
condition  with  the  aspirating  probe’s  total  pressure 
transducer  (XCQ-092).  Seven,  37,  and  111  cycle  ensemble 
averages  are  shown  superimposed  with  the  Instantaneous 
pressure  signal  where,  for  reference,  61  cycles  corresponds  to 
one  rotor  revolution.  The  seven  and  111  cycle  ensembles  are 
plotted  togethor  on  the  lower  grid  In  Figure  4-3B.  Their 
close  agreement  suggests  that  the  periodic  nature  of  the  total 
pressure  profile  can  be  found  with  relatively  few  cycles. 
Ensemble  averages  of  total  temperature  are  shown  in  Figures 
4-4  A  and  B  for  the  292  msec  test  point  with  results  similar 
to  those  seen  with  total  pressure. 

Ensemble  averages  for  total  pressure  and  temperature  are 
presented  In  the  same  format  as  above  except  at  the  360  msec 
test  point,  l.e.  during  the  traverse.  Total  pressure  ratios 
are  shown  In  Figures  4-5  A  and  B  while  the  total  temperature 
ratios  are  shown  In  Figures  4-6  A  and  B.  For  both  pressure 
and  temperature,  the  111  cycle  ensembles  are  seen  to  differ 
significantly  from  the  37  and  7  cycle  traces.  The  difference. 
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the  probe  Is  located  approximately  fn  the  core  flow  region  of 
the  upstream  vane.  In  the  second,  at  360  msec,  the  probe  Is 
located  In  the  wake  region  of  the  upstream  vane  as  suggested 
by  the  low  total  pressure  reading.  In  the  third,  at  430  msec, 
the  probe  Is  located  roughly  In  the  transition  region  between 
the  vane  core  and  wake  flow.  The  data  are  presented  In 
nond Imens lonal  form,  being  normalized  by  the  vane  Inlet 
{supply  tank)  total  temperature  and  pressure.  These  values 
change  appreciably  during  the  blowdown  run,  however,  so  the 
reference  temperature  and  pressure  must  be  corrected  for  each 
test  time.  These  corrected  reference  values  are  given  for 
each  set  of  data  presented.  Finally,  It  Is  mentioned  that  the 
data  presented  In  this  Chapter  are  the  first  published  set  of 
high-frequency  total  temperature  and  pressure  data  taken  In  a 
turbine  stage. 

4.2  Experimental  Measurements 

High-frequency  measurements  taken  In  the  Turbine  Blowdown 
Facility  by  the  aspirating  probe,  angle  probe,  and  total 
pressure  probes  are  now  presented.  First,  however,  the 
approach  used  to  ensemble  average  the  data  Is  discussed. 

4.2.1  Ensembllng  Technique 

The  data  presented  In  the  following  sections  are  ensemble 
averaged  to  reveal  the  periodic  nature  of  the  flowfleld 
structure.  For  the  turbine  design  point  corrected  speed  and 
mass  flow,  the  reduced  frequency  defined  In  Chapter  1  equals 
1.07  (Ua158  m/s,  t»0.045  m,  bx»0.025m,  and  Cx»82.3  m/s).  This 
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In  Test  1,  the  XCQ-062  total  pressure  probe  was  traversed 
circumferentially  downstream  of  the  turbine  at  the  rotor  exit 
midspan;  In  Tests  2-4,  the  aspirating  probe  was  traversed. 
The  traverse  length  was  60  mm,  or  approximately  1.33  vane 
gaps.  Both  the  probe  position  and  translation  path  for  the 
traverses  are  shown  In  Figure  4-1  along  with  a  typical  time 
trace.  The  time  trace  Indicates  the  traverse  linearity.  The 
angle  probe  was  not  traversed. 

For  later  reference,  values  for  the  turbine  midspan 
velocity  triangle,  obtained  from  a  streamline  curvature 
program,  are  given  below  In  Table  4-2.  These  apply  to  the 
100X  corrected  speed  operating  condition. 


Corrected  Speed,  NC  100X 

Blade  Speed,  wr  or  U  Cm/s]  161.9 


Rotor  Absolute  Inlet  Mach  Number,  M  a3  1.189 
Rotor  Absolute  Inlet  Angle,  0  3  73.69 
Rotor  Relative  Inlet  Mach  Number,  M  r3  1.269 
Rotor  Relative  Inlet  Angle,  B  3  59.52 
Rotor  Absolute  Exit  Mach  Number ,  M  a4  0.738 
Rotor  Absolute  Exit  Angle,  0  4  -33.88 
Rotor  Relative  Exit  Mach  Number,  M  r4  1.235 
Rotor  Relative  Exit  Angle,  B  4  -60.39 


Total  Pressure  Ratio  0.2591 

Total  Temperature  Ratio  0.7699 


Table  4-2:  Summary  of  Turbine  Midspan  Velocity 
Triangle  at  100%  Corrected  Speed 


In  the  following  section,  we  present  data  at  three  points 
during  the  tests,  as  Indicated  In  Figure  4-2,  corresponding  to 
three  traverse  positions.  In  the  first  position,  at  292  msec. 
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CHAPTER  4 

THE  RESULTS  —  FLOWFIELO  MEASUREMENTS  AND  ANALYSIS 

H Igh-f requency  measurements  taken  In  the  Turbine  Blowdown 
Facility  are  described  In  this  chapter.  After  their 
description,  the  data  are  briefly  analyzed  In  two  ways. 
First,  the  effect  of  using  t Ime-resol wed  pressure  measurements 
to  estimate  turbine  performance  (Instead  of  using  the 
ttme-averged  pressure  as  measured  by  a  Pitot  probe)  Is 
Investigated.  Second,  we  dertve  the  Euler  turbine  equation 
and  use  It  In  conjunction  with  the  angle  probe  measurements  to 
calculate  the  turbine  temperature  profile.  Subsequently,  this 
calculated  profile  is  compared  to  the  measured  temperature 
profile  taken  with  the  aspirating  probe.  Before  describing 
the  measurements  and  presenting  their  analysis,  however,  the 
test  matrix  completed  for  this  thesis  Is  outlined. 

4.1  Experiment  Summary 

The  tests  reported  In  this  thesis  are  summarized  In  Table 
4-1  below: 


Test 

Corrected 

Speed 

Initial 

Temperature,  K 

Initial 
Pressure,  atm 

1 

NC  *  100% 

475.5 

4.134 

2 

NC  -  100X 

472.8 

4.137 

3 

NC  *  125X 

486.7 

4.137 

4 

NC  -  80% 

476.1 

4.137 

Table  4-1:  Summary  of  Test  Conditions 
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3.4  Summary  of  Shock  Tube  Results 

The  shock  tube  proved  to  be  a  useful  1  and  easily 
Implemented  tool  for  determining  the  frequency  response 
characteristics  of  both  the  aspirating  probe  and  angle  probe. 
The  aspirating  probe’s  undamped  natural  frequency  was  found  to 
be  15.46  kHz  for  the  turbine  test  conditions  with  a  damping 
ratio  of  .36.  The  angle  probe’s  characteristic  frequency  was 
estimated  to  be  45.5  kHz  with  a  measured  rise  time  of  5  usee. 
These  results  confirm  each  probe’s  applicability  for  use  In 
the  turbine  facility.  Measurements  performed  by  the  probes  In 
the  turbine  facility,  along  with  their  analysis,  are  discussed 
In  the  following  chapter. 
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Is  suitable  for  use  In  the  turbine.  (It  Is  noted  In  £13, 
however,  that  this  probe's  'reduced  frequency’  Is  only  12 
kHz  )  . 

The  characteristic  frequency  above  notwithstanding,  we 
must  still  determine  the  dynamic  Interaction  of  the  probe's 
set  of  pressure  transducers.  In  practice,  the  measurements  of 
each  transducer  at  some  Instant  In  time  are  reduced  to  obtain 
flowfleld  properltles  which  are  assigned  to  that  Instant  In 
time.  Clearly,  errors  will  occur  If,  due  to  the  probe’s 
geometry  and  location  of  the  stagnation  point,  unsteady 
flowfleld  conditions  persist.  Vortex  shedding  or  Intermlttant 
flow  separation  are  examples  which  might  occur. 
Unfortunately,  no  techniques  are  currently  available  for 
generating  wel 1 -character f zed  temperature  and  pressure 
fluctuations  at  frequencies  greater  than  10  kHz  £203;  hence, 
we  cannot  Investigate  the  probes  response  to  a  fluctuating 
flowfleld.  We  can,  however,  measure  the  angle  probes  response 
to  a  step  change  In  pressure  and  In  so  doing  determine  Its 
rise  time.  These  measurements  are  now  presented. 

3.3.2  Measurement  of  Angle  Probe  Frequency  Response 

The  angle  probe  was  placed  In  the  shock  tube  and  tested 
In  air.  In  all  cases,  we  were  able  to  excite  the  natural 
frequency  of  the  angle  probe  pressure  transducers,  which  are 
on  the  order  of  200  kHz.  A  typical  pressure  transducer 
response  Is  Illustrated  In  Figure  3-8.  Rise  times  of  the 
pressure  transducers  were  found  to  be  on  the  order  of  5  usee 
with  settling  times  on  the  order  of  180  usee. 
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the  angle  probe  pressure  transducers  are  addressed  and 
followed  by  a  brief  presentation  of  the  angle  probe  rise  time 
measurements . 

3.3.1  Angle  Probe  Dynamic  Coupling 

Four  miniature  silicon  diaphragm  pressure  transducers 

(Kullte  XCQ-093)  are  mounted  on  the  angle  probe  as  shown  In 

Figure  2-7.  These  transducers  have  natural  frequencies  on  the 

order  of  200  kHz,  or  higher,  depending  on  their  rated 

pressure.  In  using  the  probe  for  flowfleld  measurements, 

however,  we  are  more  Interested  In  the  dynamic  Interaction  of 

the  entire  set  of  transducers  since  they  are  used  together  to 

retrieve  flowfleld  properties. 

A  characteristic  frequency  of  the  angle  probe  can  be 

found  based  on  the  probe  diameter  and  anticipated  flowfleld 

% 

conditions  downstream  of  the  turbine  rotor  as  shown  In  Table 
3-4  below: 


Probe  diameter,  D 

CMM3 

3.30 

Mach  Number 

0.60 

Temperature,  T 

C  K  3 

355.0 

Gamma 

1.29 

Gas  Constant,  R 

CJ/KG  K] 

137.0 

Velocity,  V 

CM/S  3 

150.3 

Character  1 st  1  c 
Frequency,  F  ■  V/D 

C kHz  3 

45.5 

Table  3-4:  Angle  Probe  Characteristic 
Frequency 


Based  on  this  simple  characteristic  frequency,  the  angle  probe 
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limiting  characteristics  associated  with  the  hot  wires 
themselves  become  Important.  The  minimum  wire 
1 ength-to-d 1 ameter  ratio  required  for  negligible  conduction 
losses  to  the  probe  body  Is  200  to  500  C  3  2  3  .  And  given  that 
the  smallest  available  hot  wires  are  2.5  to  5.0  um  In 
diameter,  the  minimum  probe  channel  diameter  becomes 
essentially  set  at  0.5  to  1.0  mm.  We  may  then  choose  the  area 
of  the  sonic  orifice,  conslstant  with  construction 
limitations,  so  as  to  give  a  low  channel  Mach  number.  Thus, 
we  can  limit  the  adverse  effect  of  the  subsonic  mean  flow  on 
the  probe  natural  frequency.  Though  the  applicability  of 
equation  3.8  Is  questionable,  It  at  least  suggests  that  a  low 
Mach  number  will  also  reduce  the  level  of  damping  Introduced 
by  the  flow. 

In  effect,  the  channel  length  of  the  probe  Is  the  primary 
variable  controlling  the  natural  frequency.  Consideration 
must  be  given  to  the  Inlet  design  to  ensure  good  off-angle 
performance  and  a  smooth  diffusion  of  the  flow  Into  the 
constant  area  channel.  But  In  theory,  It  seems  the  designer 
need  only  choose  a  length  commensurate  with  construction 
constraints  to  obtain  the  desired  natural  frequency.  In  the 
present  case,  where  we  would  like  a  natural  frequency  of  40  to 
60  kHz,  we  would  need  to  reduce  the  probe’s  channel  length  by 
one-half  or  more. 

3.3  Angle  Probe  Frequency  Response 

Frequency  response  aspects  of  the  angle  probe  are 
discussed  In  this  section.  In  particular,  dynamic  coupling  of 
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unreasonable,  though,  because  at  the  blade  passing  frequency 
It  yields  a  magnification  ratio  of  only  1.2,  approximately. 

In  summary,  the  natural  frequency  and  damping  ratio  of 
the  aspirating  probe  have  been  measured  In  both  air  and 
freon-12.  The  results  confirm  the  applicability  of  equation 
3.7  for  predicting  undamped  natural  frequency,  but  suggest 
that  equation  3.8,  the  damping  ratio  prediction,  Is  not 
applicable.  From  these  results,  the  aspirating  probe’s 
predicted  characteristics  for  the  argon-f reon- 1 2  mixture  used 
In  the  TBF  are  f*15.36  kHz  and  zb0.36.  The  damping  ratio  Is 
an  approximation  based  on  the  measured  results.  It  Is  felt 
that  these  values  of  natural  frequency  and  damping  ratio  are 
sufficient  to  measure  the  b 1 ade-to-b 1 ade  temperature  and 
pressure  distributions  downstream  of  a  turbine  rotor.  Given 
these  results,  methods  for  Improving  the  aspirating  probe 
frequency  response  are  now  considered. 

3.2.4  Design  Considerations  for  Frequency  Response  Improvement 

The  aspirating  probe  natural  frequency  Is  a  function  of 
Its  channel  length,  the  sonic  speed,  and  the  channel  Mach 
number,  as  was  shown  In  equation  3.7.  The  sonic  speed  Is  a 
function  of  the  test  gas  and  temperature,  and  therefore,  more 
an  application  than  design  consideration.  So  on  a  next 
Iteration  of  the  aspirating  probe  design,  we  are  left  with 
only  the  channel  length  and  channel  Mach  number  (a  function  of 
channel-to  on  1 c-or If  1 ce  area  ratio)  as  design  variables  for 
Improving  the  natural  frequency. 

In  considering  a  design  for  Improved  frequency  response. 


PAGE  69 


0.32,  ts  shown  In  Figure  3-7. 

We  may  now  predict  the  aspirating  probe’s  frequency 
response  for  the  TBF’s  argon-freon- 1 2  test  gas  mixture. 
Results  for  a  freon-12  mass  fraction  of  0.511%  are  shown  In 
Table  3-3  below: 


Test  Gas 

Argon 

-  Freon-12 

Channel  length,  L 

CMMI 

3.886 

Temperature,  T 

CK] 

355.0 

Freon  mass  fraction 

[%J 

0.511 

Gas  constant,  R  CJ/KG  K1 

137.0 

Gamma 

1 .29 

Sonic  speed,  C 

CM/S] 

250.4 

Channel  Mach  number. 

M 

0.2 

Natural  Frequency,  f 

CkHz] 

15.46 

Damping  Ratio,  z 

0.36 

(approximated  from 

measurements ) 

Table  3-3:  Aspirating  Probe  Natural  Frequency 
and  Damping  Ratio  Prediction  In  an 
Argon-Freon- 1 2  Mixture 


The  prediction  above  Is  based  on  the  typical  DC  value  for 
temperature  downstream  of  the  turbine  rotor.  The  damping 
ratio  shown  Is  an  approximation  based  on  the  values  measured 
for  air  and  freon-12.  The  natural  frequency  of  15.46  kHz  Is 
well  above  blade  passing  frequency,  but  ts  not  nearly  high 
enough  to  Investigate  wake  structure.  It  Is  felt,  though, 
that  this  natural  frequency  and  damping  ratio  pair  are 
sufficient  to  dlstlnglush  b 1 ade- to-b 1 ade  variations  In 
temperature  and  pressure.  The  damping  ratio  of  0.36  Is  rather 
less  than  the  optimum  value  of  0.6  to  0.7.  It  Is  not 


The  predicted  and  measured  undamped  frequencies  agree 
very  well,  thus  giving  credibility  to  equation  3.7.  The 
variance  In  the  freon-12  measurement  and  prediction  Is  most 
likely  associated  with  the  presence  of  air  In  the  driven  tube 
during  the  test.  The  air  would  alter  the  driven  gas 
properties,  resulting  In  a  higher  natural  frequency  similar  to 
that  measured.  Unf or tuneate 1 y ,  the  damping  ratios  do  not 
agree.  The  prediction  gives  a  damping  ratio  close  to  the 
critical  value  of  2=1;  a  value  for  which  we  would  expect  to 
see  no  oscillations  In  the  signal.  Clearly  the  damping  ratio 
prediction,  equation  3.8,  Is  not  applicable  to  the  aspirating 
probe . 

The  frequency  response  of  a  simple  hot  wire  similar  to 
that  used  In  the  asplratfng  probe  was  also  measured  In  the 
shock  tube.  Tests  were  accomplished  In  air  for  the  same 
diaphragm  pressure  ratio  and  temperature  as  used  with  the 
aspirating  probe.  This  was  done  to  ensure  the  measured 
frequencies  were  a  function  of  the  aspirating  probe  geometry 
and  gas  properties  alone,  and  not  associated  with  the  hot  wire 
circuitry  or  shock  tube  Itself.  Based  on  the  stability 
setting  of  the  anemometer  circuit,  damped  natural  frequencies 
of  45.0  to  65.0  kHz  were  observed  for  the  simple  hot  wire. 
These  values  are  well  above  those  measured  for  the  aspirating 
probe  and  further  tend  to  verify  Whiteheads  result  for 
predicting  natural  frequency.  A  typical  hot  wire  response,  a 
natural  undamped  frequency  of  59  kHz  and  damping  ratio  of 
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frequency  Is  always  less  than  the  undamped  natural  frequency. 
We  can,  however,  extract  the  required  Information  using  the 
following  relationships  and  definitions  given  In  Figure  3-6 
C  3 1  ]  > 


[3.101 


[3.111 


Using  the  equations  above  and  the  measurements  In  Figures  3-4 
and  3-5  ,  we  obtain  the  results  shown  In  Table  3-2  below. 


Test  Gas 

A 1  r 

Freon- 1 2 

a 

[Volts] 

0.52 

1 .28 

A 

t Vo  1 ts 1 

0.  14 

0.44 

Period,  T 

[Ml crosec 1 

55.0 

95.0 

Pr ed 1 cted 

damp  1 ng 

ratio,  z , 

from  3.8 

0.97 

0.89 

Measured 

damp  1 ng 

rat 1o,  z , 

from  3.10 

0.39 

0.32 

Pred 1 cted 

undamped 

natural  frequency,  f 

from  3.7 

[kHz! 

21.63 

9.53 

Measured 

undamped 

natural  frequency,  f 

from  3.11 

[kHz  1 

21.35 

11.11 

Table  3-2»  Comparison  of  Predicted  to  Measured 
Aspirating  Probe  Natural  Frequency 
and  Damping  Ratio 
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With  the  natural  frequency  and  damping  ratio  predicted, 
the  task  Is  now  to  measure  these  characteristics  of  the  probe 
and  compare  the  measurements  to  the  predictions.  If  they 
agree,  we  may  use  equations  3.7  and  3.8  to  predict  the  probe’s 
frequency  response  In  argon-f reon- 1 2 . 

The  aspirating  probe  was  mounted  In  the  shock  tube  test 
facility  and  subjected  to  step-function  changes  In  both 
pressure  and  temperature.  Separate  tests  were  accomplished  In 
both  air  and  freon-12  at  diaphragm  pressure  ratios  of 
approximately  1.5,  yielding  pressure  ratios  at  the  probe  of 
about  1.25.  This  level  of  pressure  change  Is  similar  to  what 
the  probe  would  experience  downstream  of  the  turbine  rotor  In 
the  TBF .  Typical  probe  responses  to  the  step  change  In 
properties,  along  with  their  respective  power  spectrums,  are 
given  In  Figures  3-4  A  and  B  for  air  and  Figures  3-5  A  and  B 
for  freon-12.  Test  conditions  are  noted  on  each  figure. 

We  can  determine  the  aspirating  probe’s  natural  frequency 
and  damping  ratio  from  Figures  3-4  and  3-5.  First,  It  must  be 
realized  that  a  systems  undamped  natural  frequency  cannot  be 
observed  experimentally  unless  the  system  has  zero  damping. 
What  Is  actually  observed  Is  the  damped  natural  frequency 
which  Is  related  to  the  undamped  natural  frequency  as  follows 
C  30]  : 


w  damped 


w 


natural 


[3.9] 


This  result  applies  only  to  underdamped  systems  where  z  <  1.0. 
From  equation  3.9  It  Is  clear  that  the  experimentally  observed 
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temperature.  No  attempt  was  made  to  determine  the  flowfleld 
total  temperature  and  pressure  with  the  hot  wires  alone.  On 
each  figure  listed  In  Table  4-3,  the  total  pressure  ratio  and 
total  temperature  ratio  measured  from  each  wire  are  plotted 
against  time.  For  the  100%  corrected  speed  condition,  an 
additional  figure  Is  given  which  compares  the  ensembled  and 
instantaneous  temperature  signals.  The  second  time  point  (360 
msec)  for  each  speed  condition  varies  slightly.  This  Is  due 
to  differences  In  the  probe  traversal  from  test  to  test--the 
probe  position  is  the  same  for  each  data  set. 

In  theory,  the  temperature  reduced  from  each  hot  wire 
should  be  the  same.  In  practice,  they  vary  slightly.  Maximum 
differences  in  temperature  between  the  two  wires  Is  on  the 
order  of  18  C,  but  these  occur  In  the  peaks  and  valleys  of  the 
ensemble  averages.  Generally,  the  form  of  the  ensembled 
temperature  traces  agree  very  well.  The  temperature 
sensitivities  of  the  wires  are  different,  as  mentioned  In 
Chapter  2,  but  the  calibration  data  should  negate  this  effect. 
Except  that  the  wires  are  spatially  separated  (radially  In 
reference  to  the  blade  trailing  edge,  and  only  1%  In  blade 
span),  no  explanations  for  the  differences  In  probe  output  are 
given  here. 

For  the  100%  corrected  speed  condition,  the  temperature 
ratio  Is  seen  to  have  a  characteristic  frequency  twice  that  of 
blade  passing.  This  is  most  clearly  evidenced  (n  the  292  msec 
and  430  msec  test  points.  The  peak-to-peak  variation  Is 
typically  20%,  or  65  C.  The  ensembles  from  these  test  points 
compare  very  favorably  with  their  Instantaneous  temperature 
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traces.  During  the  traverse,  however,  the  temperatures  are 
not  so  cleanly  structured.  Ensembling  the  data  over  seven 
cycles  Instead  of  37  produces  the  same  result,  as  is  shown  in 
Figure  4-7B.  This  suggests  there  is  not  a  problem  with  the 
ensembling  technique;  rather,  the  temperature  in  the  vane 
wake  region  has  a  more  random  structure.  This  random  nature 
is  seen  in  Figure  4-10,  where  it  is  clear  that  the  ensemble 
averages  do  not  agree  well  with  the  Instantaneous 
measurements.  For  all  the  test  points,  the  peak  in  total 
temperature  ratio  corresponds  roughly  with  the  trough  in  total 
pressure  ratio.  And  assuming  the  wake  is  signalled  by  low 
total  pressure,  the  data  suggest  the  wakes  are  hot. 

For  the  125X  operating  condition,  we  see  again  the  twice 
blade-passing  frequency  of  the  temperature  profile.  The 
amplitude  of  the  second  peak  is  less  pronounced,  however,  than 
in  the  100X  speed  case.  The  maximum  peak-to-peak  temperature 
variation  for  this  speed  is  roughly  25X,  or  70  C.  The  phase 
relationship  between  the  pressure  and  temperature  has  changed 
from  the  100X  speed  condition.  Now,  the  trough  in  total 
temperature  lines  up  more  closely  with  the  trough  in  total 
pressure . 

Finally,  for  the  80X  operating  condition,  we  see 
temperature  profiles  which  are  much  less  conslstant  than  those 
from  the  other  operating  conditions.  The  profiles  contain 
frequency  components  at  three  times  blade-passing  and  greater 
with  a  different  structure  at  each  test  point.  Further 
correlation  is  required  for  this  data  set,  but  it  has  not  been 
accomplished  for  this  thesis.  Typical  peak-to-peak 
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temperature  variation  is  11%,  or  37  C.  And  like  the  100% 
corrected  speed  case,  the  temperature  ratio  is  seen  to  vary 
alternately  with  the  pressure  ratio. 


4.2.3  Angle  Probe  Measurements 

Total  and  static  pressure  ratios,  Mach  numbers,  and 
tangential  flow  angles  reduced  from  the  angle  probe  pressure 
measurements  are  presented  In  this  section.  These 
measurements  are  used  later  with  the  Euler  turbine  equation  to 
predict  the  rotor  exit  temperature  profile.  Unf ortuneatel y , 
during  the  course  of  testing,  the  angle  probe’s  transducer 
number  4  became  1 nter m 1 ttant .  As  a  result,  no  radial  flow 
Information  was  obtained.  The  angle  probe  was  not  traversed 
during  the  tests.  Its  stationary  position  corresponds  to  the 
360  msec  test  point  for  the  traversed  probes.  Finally,  only 
the  measurements  from  the  100%  corrected  speed  condition  are 
presented  here. 

Both  the  absolute  and  relative  tangential  flow  angles  are 
plotted  In  Figure  4-19.  The  relative  angle  Is  calculated 
using  the  absolute  Mach  number  and  tangential  angle  measured 
by  the  angle  probe,  along  with  the  wheel  speed.  The  probe  was 
offset  from  axial  In  the  laboratory  frame  by  34  degrees  so  as 
to  face  directly  Into  the  flowfleld.  The  streamline  curvature 
program  mentioned  In  Section  4-1  predicts  values  of  34.88  and 
-60.39  for  the  absolute  and  relative  tangential  angles 
respectively.  From  the  Figure,  It  Is  seen  that  the  absolute 
tangential  flow  angle  has  a  peak-to-peak  variation  of  30 
degrees,  with  the  mean  value  corresponding  to  the  streamline 
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curvature  prediction.  Likewise,  the  relative  tangential  angle 
has  a  significant  peak-to-peak  variation  of  15  degrees.  This 
Is  Interesting  since  it  Is  normally  assumed  that  the  relative 
tangetial  flow  angle  is  constant  except  for  a  small  deviation 
angle  at  the  blade  trailing  edge. 

The  total  pressure  ratio  and  absolute  Mach  number  are 
shown  in  Figure  4-2 0.  In  Figures  4-21  and  4-22,  the  axial 
Mach  number  and  tangential  Mach  number,  respectively,  are 
shown  with  the  total  pressure  ratio.  And  In  Figure  4-23,  the 
axial  and  tangential  Mach  number  components  are  plotted 
togethor .  From  Table  4-2,  the  predicted  absolute  exit  Mach 
number  Is  0.738  for  the  blade  midspan.  This  Indicates  the  DC 
value  for  the  measured  absolute  Mach  number  Is  low  by 
approximately  10%.  From  Figure  4-20,  we  see  the  trough  in 
absolute  Mach  number  lines  up  with  the  minimum  total  pressure 
as  one  would  expect.  The  Mach  number  trace,  though  generally 
increasing  with  time,  .l.e  from  the  blade  pressure  to  suction 
side,  Is  very  non-uniform.  This  aspect  of  the  Mach  number 
profile  was  unexpected  and  requires  further  analysis.  The 
axial  and  tangential  Mach  number  components  are  seen  to  be  180 
degrees  out  of  phase,  with  the  minimum  In  axial  Mach  number 
occurring  at  the  minimum  total  pressure.  In  contrast,  the 
maximum  tangential  Mach  occurs  at  the  pressure  minimum. 

Finally,  the  total  pressure  ratio  and  static  pressure 
ratio  are  shown  together  In  Figure  4-24.  Surprisingly,  the 
static  pressure  has  roughly  the  same  blade-  to-blade  profile 
as  the  total  pressure.  Clearly,  the  static  pressure  Is  not 
constant  through  the  wake  region  as  Is  normally  assumed. 
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4.2.4  Total  Pressure  Probe  Measurements 

The  Kullte  XCQ-062  is  the  smallest  pressure  transducer 
available  which  provides  accurate  DC  and  AC  pressure 
measurements.  A  total  pressure  probe  containing  this 
transducer  yields  the  pressure  measurements  shown  In  Figure 
4-25.  In  this  figure,  both  the  ensembled  and  Instantaneous 
traces  are  shown.  On  the  figure’s  lower  grid,  the  ensembled 
pressure  ratio  measured  by  the  probe  is  superimposed  on  an 
ensembled  trace  taken  by  the  aspirating  probe’s  XCQ-093 
transducer.  The  measurements  were  taken  In  the  same  position 
and  test  time  but  in  different  tests.  The  form  of  the 
ensembled  pressure  traces  are  quite  similar.  Each  show  a 
strong  pressure  fluctuation  at  the  blade  passing  frequency 
with  a  smaller  amplitude  frequency  superimposed.  This 
secondary  frequency  is  most  likely  associated  with  the  blade 
passage  vortex.  The  time-resolved  pressure  measurement 
obtained  with  the  XCQ-062  is  employed  in  the  following  section 
to  investigate  turbine  performance. 

4.3  Data  Analysis 

In  this  section,  the  practice  of  using  time  averaged 
pressures  to  calculate  turbine  performance  Is  discussed  in 
light  of  that  component’s  characteristically  unsteady 
flowfield.  Comparison  Is  made  between  the  time-average  of  the 
time-resolved  pressure  signal  described  above,  and  the 
time-averaged  pressure  signal  as  might  be  measured  by  a  Pitot 
probe.  Next,  the  nond  imens (ona 1 ized  Euler  turbine  equation  is 
derived  and  used  to  show  the  dependence  of  the  stage 
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temperature  ratio  on  the  blade  turning  angles,  flow  Mach 
numbers  and  rotor  speed.  Finally,  comparison  Is  made  between 
the  temperature  profile  calculated  from  the  Euler  turbine 
equation  and  that  measured  by  the  aspirating  probe. 

4.3.1  Fluctuating  Pressure  Influence  on  Performance 

Calculations 

Adiabatic  efficiency  serves  as  a  measure  of  turbine 
performance  and  can  be  defined  as  follows* 


where  Tt  Is  the  gas  total  temperature,  Pt  Is  the  total 
pressure,  and  3  and  4  refer  to  the  rotor  Inlet  and  exit 
respect  1 vel y .  This  formulation  applies  only  for  a  mass-flow 
weighted  average  of  the  properties  as  indicated  by  the  carats. 
If  the  variations  In  flow  velocity  and  density  are  small  then 
a  time-weighted  average  of  the  properties  may  be  used  as  an 
approximation  to  the  mass-averaged  values.  Time-weighted 
averages  will  be  used  here. 

Together,  a  Pitot  tube,  connecting  lines,  and  an  Inert 
measuring  Instrument  comprise  an  ’Integrating  pneumatic 
system’  which  can  measure  time-averaged  pressures  In 
fluctuating  flows  such  as  those  downstream  of  a  turbine  rotor. 
In  practice,  such  systems  are  commonly  used  In  turbine  test 
facilities  to  measure  the  time-averaged  pressure  which  Is  then 
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used  with  Equation  4.1  to  estimate  the  turbine  performance. 
However,  as  shown  by  Weyer  133],  measurements  of  this  sort  are 
subject  to  large  errors  which  depend  upon  the  strength  of  the 
pressure  oscillations,  the  pressure  wave  shape,  and  probe 
construction.  Figure  4-26  shows  the  devfatlon  of  an  ordinary 
Pitot  probe  reading  from  the  correct  time-  averaged  pressure 
measured  downstream  of  a  transonic  axial  flow  compressor.  In 
this  case,  errors  of  up  to  8%  occur  as  the  rotor  speed,  and 
hence  flowfleld  fluctuations.  Increase  (temperature 
measurement  errors  have  also  been  noted  but  these  are  not 
addressed  here).  If  taken  on  face  value,  the  erroneous 
measurements  would  suggest  a  higher  efficiency  than  the 
compressor  actually  possesses.  We  will  Investigate  here  the 
magnitude  of  similar  errors  which  can  occur  in  turbines. 

The  measuring  system  model  described  above,  one 
consisting  of  a  throttle  point,  connecting  duct,  and  Inert 
measuring  Instrument,  will  sense  some  time  averaged  pressure, 
Pbar ,  when  Immersed  In  a  flowfleld  with  fluctuating  pressure, 
Pl(t).  (Basically,  Pbar  Is  the  ’probe  averaged’  pressure).  A 
schematic  of  this  measuring  system  Is  shown  In  Figure  4-27.  A 
relationship  between  the  mean  pressure  sensed  by  the  measuring 
system,  Pbar,  and  the  actual  time-varying  pressure,  Pl(t),  Is 
given  by  Weyert 
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where  2n  Is  the  number  of  Identical  time  steps  In  a  sampled 
time  Interval,  and  k  Is  the  number  of  samples  for  which  Pitt) 
Is  greater  than  Pbar  .  As  suggested  by  the  definition  of  k, 
the  first  term  In  Equation  4.2  applies  when  Pl(t)  Is  greater 
than  Pbar,  the  second  when  Pitt)  Is  less  than  Pbar.  Weyer * s 
derivation  of  this  relationship  Is  based  essentially  upon  the 
application  of  the  Mass  Conservation  Principle  across  the 
throttle  plate  shown  in  Figure  4-27,  and  Is  fully  detailed  In 
C  34  3  . 

Having  measured  the  unsteady  pressure  downstream  of  the 
test  turbine,  we  can  now  time  average  It  directly.  Next,  we 
can  use  the  same  time-resolved  signal  to  calculate  Pbar  from 
Equation  4.2,  with  the  result  shown  In  Figure  4-28  1351.  For 
this  calculation,  the  pressure  signal  from  the  XCQ-062 
pressure  transducer  Is  employed  with  averaging  done  over  each 
blade  passing  cycle.  The  difference  between  Pbar  and  the 
correct  time-averaged  signal,  Pavg,  Is  also  shown  In  Figure 
4-28.  Their  difference  Is  seen  to  fluctuate  between  _+  2. AX, 
with  a  mean  difference  of  ^0.4X. 

A  streamline  curvature  program  gives  for  the  turbine 
temperature  and  pressure  ratios  0.7699  and  0.2591, 
respectively,  at  100X  corrected  speed.  With  these  ratios  and 
1.3  for  the  test  gas  ratio  of  specific  heats,  we  can  use 
Equation  4.1  to  calculate  the  turbine  adiabatic  efficiency. 
The  result  Is  86. 0X.  Using  the  maximum  difference  between 
Pbar  and  Pavg,  we  see  that  errors  of  1.25  points  In  adiabatic 
efficiency  are  possible.  Realistically,  though,  the  error 
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associated  with  the  mean  pressure  difference  Is  more  likely  to 
be  realized.  This  error  Is  only  0.25  points  In  adiabatic 
efficiency.  It  seems  small,  but  Is  actually  significant  given 
the  high  premium  on  turbine  adiabatic  efficiency. 

4.3.2  Euler  Turbine  Equation 

The  Euler  turbine  equation  describes  the  energy  transfer 
process  In  turbomachines  and  Is  obtained  by  equating  the 
energy  and  momentum  equations  via  the  turbine  power  and 
torque.  For  the  case  of  steady,  two-dimensional,  axlsymmetrlc 
flow  the  Euler  turbine  equation  takes  the  following  form  with 
the  velocity  triangle  paramters  defined  In  Table  4-2: 


cp4  T 4  -  cp3  T3  -  w  <r4  Vo4  -  r3  Vo3 )  14.31 


where  cp  is  the  gas  specific  heat  (a  function  of  temperature), 
T  the  gas  total  temperature,  w  the  rotor  rotational  speed,  r 
the  streamtube  radius  (which  may  shift  with  passage  through 
the  blade  row),  and  Vo  the  tangential  flow  velocity.  Stations 
upstream  and  downstream  of  the  turbine  rotor  are  denoted  by  3 
and  4  respectively. 

We  can  redefine  the  rotor  speeds,  wr3  and  wr4,  and 
tangential  flow  velocities,  Vo3  and  Vo4,  In  Equation  4.3  In 
terms  of  Mach  number: 
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where  M  fs  the  Mach  number  based  on  the  local  gas  temperature, 
R  the  gas  constant,  and  v  the  ratio  of  specific  heats.  The 
blade  tip  Mach  number,  Mtlp,  is  based  on  the  vane  inlet 
conditions.  Incorporating  these  into  Equation  4.3  yields: 
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For  a  given  set  of  flowfield  conditions,  the  first  term 
in  the  right  hand  side  of  Equation  4.7  can  be  considered  a 
constant;  the  second  a  constant  times  the  square  root  of  the 
temperature  ratio.  After  some  manipulation,  the  equation  can 
be  reduced  to  the  quadratic  form  from  which  the  temperature 
ratio,  T4/T3,  is  readily  calculated.  Clearly,  we  see  the 
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turbine  temperature  ratfo  Is  a  function  of  rotor  speed,  Inlet 
and  exit  flow  angle,  and  Inlet  and  exit  Mach  number. 

We  can  predict  the  rotor  exit  temperature  profile  with 
Equation  4.7  by  using  the  Mach  number  and  tangential  flow 
angle  measured  by  the  angle  probe.  Unf or tuneate 1 y ,  no  vane 
exit  measurements  were  taken,  so  the  rotor  Inlet  Mach  number 
and  tangential  angle  must  be  determined  and  Input.  These 
values,  given  In  Table  4-2,  are  obtained  with  a  streamline 
curvature  program.  Input  directly,  and  held  constant. 

The  result  of  this  calculation,  the  Euler  total 
temperature  ratio.  Is  shown  In  Figure  4.29  along  with  the 
total  pressure  ratio  measured  by  the  angle  probe.  The  Euler 
temperature  profile  varies  at  the  blade  passing  frequency, 
unlike  the  measured  profile  which  varies  at  twice  that 
frequency.  The  amplitude  of  each  Is  roughly  the  same,  but  the 
Euler  profile  Is  DC  shifted  upward  20%  from  the  measured 
values.  Finally,  and  again  unlike  the  measured  temperature, 
the  Euler  temperature  ratio  varies  roughly  In  phase  with  the 
pressure  ratio,  with  the  lowest  total  temperature  occurring 
approximately  with  the  lowest  total  pressure.  These 
discrepancies  are  only  noted  and  not  analyzed  further  In  this 
thesis.  Their  occurrence,  however.  Is  not  surprising  since 
the  Euler  turbine  equation  as  derived  only  applies  to  steady, 
two-dimensional,  axlsymmetrlc  flow. 
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CHAPTER  5 
THE  CONCLUSIONS 

Time-resolved  flowfleld  measurements  have  been 
successfully  taken  In  the  MIT  Turbine  Blowdown  Facility.  The 
properties  measured  Include  total  temperature,  total  pressure, 
static  pressure,  Mach  number,  and  tangential  flow  angle.  The 
original  objectives  were:  to  measure  the  b 1 ade-to-b 1 ade  total 
temperature  profile  for  subsequent  compar Iston  with  a 
prediction  by  the  Euler  turbine  equation;  to  determine  the 
effect  of  using  time-averaged  pressures  to  calculate  turbine 
performance;  and  to  provide  a  complete  set  of  time-resolved 
turbine  stage  data.  No  time-resolved  turbine  data  are 
currently  available  In  the  literature.  A  preliminary 
objective  was  to  determine  the  frequency  response  of  the 
Instrumentation  employed  to  make  the  flowfleld  measurements. 
With  these  objectives  In  mind,  the  major  results  and 
conclusions  of  this  thesis  are  now  summarized. 

First,  we  found  the  shock  tube  to  be  an  easily 
Implemented  tool  fo.  1 nvest 1 ga t I ng  Instrumentation  frequency 
response  characteristics.  In  particular,  for  the  TBF  test 
gas,  the  aspirating  probe  was  found  to  have  a  natural 
frequency  of  15.5  kHz  with  a  damping  ratio  of  0.35.  The  angle 
probe  was  found  to  have  a  rise  time  of  5  ysec  with  a  settling 
time  of  18  usee.  These  character  1  st  1  cs  of  the  aspirating 
probe  and  angle  probe  were  satisfactory  for  attempting 
flowfleld  measurements  In  the  turbine  facility. 

Second,  the  TBF  and  Its  complement  of  high-frequency 
response  Instrumentation  allowed  for  detailed  flowfleld 
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measurements  in  a  rigorously  scaled  test  environment.  As 
mentioned  above,  time-resolved  measurements  of  total 
temperature,  total  pressure,  static  pressure,  Mach  number,  and 
tangential  flow  angle  were  successfully  obtained. 

Third,  the  time-resolved  pressure  measured  with  a  Kullte 
XCQ-062  pressure  transducer  was  employed  to  investigate 
turbine  performance  calculations.  Specifically,  errors 
associated  with  using  the  time-averaged  pressure  from  a 
standard  pitot  probe,  Instead  of  the  true  time-average,  were 
shown  to  be  on  the  order  of  0.25  to  1.25  points  in  adiabatic 
ef f i c 1 ency . 

Finally,  the  turbine  exit  total  temperature  was 
calculated  using  the  Euler  turbine  equation  in  conjunction 
with  the  tangential  angle  and  Mach  number  from  the  angle 
probe.  The  predicted  profile  varied  only  at  the  blade  passing 
frequency,  unlike  the  measured  profile  which  had  a 
characteristic  frequency  of  twice  blade  passing.  Further,  the 
predicted  profile  was  DC  shifted  upward  10X,  although  its 
peak-to-peak  temperature  change  was  roughly  the  same  as  that 
measured . 
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Figure  2-4:  Schematic  of  Aspirating  Probe/Total  Pressure 
Probe  Assembly;  Cross-Section  of  Aspirating 
Probe 
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Figure  2-5At  Prediction  of  Aspirating  Probe  Calibration 
Data  With  Uncorrected  Fluid  Properties. 
Overheat  Ratio  ■  1.8 
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F 1 gure  2-5C  : 


Prediction  of  Aspirating  Probe  Calibration 
Data  With  Corrected  Fluid  Properties. 
Overheat  Ratio  ■  1.8 
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Figure  2-50:  Prediction  of  Aspirating  Probe  Calibration 
Data  With  Corrected  Fluid  Properties, 
Overheat  Ratio  =  2.B 
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Figure  2-6:  Aspirating  Probe  Combined  Calibration 
Map  ( from  [ 1 0  1  > 
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Figure  2-10D:  Angle  Probe  Coefflcent  Prediction, 

CPI  Versus  Mach  Number  with  Phi  ■  0.. 
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Figure  2-10C:  Angle  Probe  Coefflcent  Prediction, 

KP2  Versus  Mach  Number  with  Phi  *  0.0 
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DETERMINE  Pg  FROM  Kp2  AND/OR  Kp3 


DETERMINE  PT  FROM  Cp|  and  H23 

i 

i 

DO  PT  4  Pg  CONVERGE? 

J/yes 

CALCULATE  M,  Mx,  M^,  M0 


i  ~ 

DOES  M  CONVERGE? 

I  YES 

Pj,  Pg,  e,  M,  M0,  Mx 


e  2-12:  Flowchart  of  Angle  Probe  Data 
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Figure  4-6Bi  Ensemble  Average  of  Total  Temperature  Ratio  with 
111  Cycle,  NC*100X;  Comparison  of  111  and  7 
Cycle  Ensemble  Averages,  Tlme'360  msec 
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Figure  4-5Bi  Ensemble  Average  of  Total  Pressure  Ratio  with 
111  Cycle.  NC-lffBX:  Comparison  of  37  and  7 
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Figure  4-5At  Ensemble  Average  of  Total  Pressure  Ratio  with 
7  and  37  Cycles:  NC»100X,  T1me=360  msec 
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Figure  4-*B:  Ensemble  Average  of  Total  Temperature  Ratio  with 
111  Cycle.  NC»100%;  Comparison  of  111  and  7 
Cycle  Ensemble  Averages.  Tfme*292  msec 
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Figure  4-4A;  Ensemble  Average  of  Total  Temperature  Ratio  with 
7  and  37  Cycles;  NC»100%,  Tlme«29Z  msec 
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Figure  4-3Bi  Ensemble  Average  of  Total  Pressure  Ratio  with 
111  Cycle,  NC=100%;  Comparison  of  111  and  7 
Cycle  Ensemble  Averages,  Ttme=Z92  msec 
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Figure  4  —  2 1 t  Total  Pressure  Ratio  and  Axial  Mach  Number 
Versus  Time 
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Figure  4-25:  Comparison  of  Total  Pressure  Ratio  as 

Measured  by  XCQ-093  and  XCQ-063  Silicon 
Diaphragm  Pressure  Transducers 
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